I 


w w i 


/ 


AFFDL-TR-73-42 

"■ 20070917040 


DEVELOPMENT  AND  EVALUATION  OF  METHODS  OF 
PLANE  STRESS  FRACTURE  ANALYSIS 
Part  III.  Application  of  the  Residual  Strength  Prediction 
Technique  to  Complex  Aircraft  Structure 


NOR  THROP  CORPORA  TION 
AJRCRA  FT  DIVISION 


AUGUST  1977 


TECHNICAL  REPORT  AFFDL-TR-73-42.  PART  III, 
Report  for  Period  January  1975  - June  1975 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


Best  Available  Copy 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  Incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Govern- 
ment may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto . 

This  report  has  been  reviewed  by  the  Office  of  Information  (01)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) . At 
NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


G.P.  SENDECKYJ 
Project  Engineer 


ROBERT  M.  BADER,  Chief 

Structural  Integrity  Br 
Structural  Mechanics  Division 


Chief,  Structural  Mechanics  Division 


Conies  of  this  report  should  not  be  returned  unless  return  is  required 
by  security  considerations,  contractual  obligations,  or  notice  on  a specific 


Affi  FORCE/56760/17  October  1977  - 500 


UNCLASSIFIED, 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWhtn  Omit  Bn tnrtd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 REPORT  number 

AFFDL-TR- 73-42,  Part  111 

2,  GOVT  ACCESSION  NO* 

9.  recipient's  catalog  number 

*.  title  (and  Submit)  DEVELOPMENT  AND  EVALUATION  OF 
METHODS  OF  PLANE  STRESS  FRACTURE  ANALYSIS  Part  III 
Application  of  the  Residual  Strength  Prediction 
Technique  to  Complex  Aircraft  Structure 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Part  III  (Phase  III)  Final 
Report  - -Jan , 1975  - June  1975 

«.  PERFORMING  ORG.  REPORT  NUMBER 

NOR  75-86 

7.  author*-#) 

M * M , Ratwani  and  D . P , W 1 Ihem 

«.  CONTRACT  OR  GRANT  NUMBER!*; 

F33615-72-C-1796 

s.  performing  organization  name  and  address 

Northrop  Corporation,  Aircraft  Division 

3901  West  Broadway,  Hawthorne,  California  90250 

10.  program  element  project,  task 

AREA  4 WORK  UNIT  NUMBERS 

Project  486U 

Task  0204 

1 i.  controlling  office  n ame  and  address 

Air  Force  Flight  Dynamics  Laboratory 

W right -Patterson  Air  Force  Base,  Ohio  45433 

12.  report  date 

August  1977 

1).  NUMBER  OF  PAGES 

148 

14,  MONITORING  AGENCY  NAME  ft  ADDRESSfJf  different  from  Controlling  Office) 

if.  SECURITY  CLASS,  (of  thle  report ) 

Unc lassi f ied 

IS#,  DECLASSIFICATION/ DOWN  GRADING 

SCHEDULE 

it,  01  ST  R!  BU  T|  ON  STATEMENT  (at  thle  Report) 

Approved  for  Public  Release;  Distribution  Unlimited* 

17,  DISTRIBUTION  STATEMENT  (at  the  mbetreci  ant  trad  In  Block  70*  It  different  ham  Report) 

it  supplementary  notes 

19,  KEY  WORDS  (Continue  on  r$vene  tt  neceemery  end  identify  by  block  number) 

Structural  Residual  Strength  P rand 1 1 - Reus s Material  Behavior  Fini  te  Element 

Fracture  Mechanics  Analysis  J -Integral  Modeling 

Elastic -P last ic  Analysis  /Jr  Resistance  Curve 

Dugdale  Model  Plane  Stress  Fracture 

Biaxial  Loading  Mixed  Mode  Fracture 

20  ABSTRACT  (Continue  on  rev* tee  mid*  it  necefery  *nd  Identity  by  bJk»eA  numb*/') 

Using  the  residual  strength  technique  presented  in  Part  1L  of  this  report, 
uniaxial ly  loaded,  angle  stiffened  wing  panels  were  analyzed  and  fracture 
strengths  determined*  Excellent  correlation  was  obtained  between  experimental 
and  analytical  data  for  an  Initially  intact  and  broken  central  stringer  for  a 
six  bay  aluminum  panel  with  thin  skin,  a thick  (0.193  inch)  skin  aluminum 
panel,  and  an  all  titanium  panel*  Both  crack  arrest  and  fracture  could  be 

(See  over) 

DD  H73  EDITION  OF  l NOV  6S  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (When  Bmtm  Enlmrtd) 


UNCLASS IF IED 

SECURITY  CLASSIFICATION  OF  THIS  PAOE(Whm  Dmtm  Enfr*4) 


20  * ABSTRACT  (Continued) 

predicted  using  the  tangency  conditions  between  the  crack  driving  force  curves 
(/J)  and  the  crack  growth  resistance  curve  i/JrJ,  The  influence  of  biaxial 
load  ratio  (tension  fields)  on  crack  opening  displacement,  plastic  zone  size  and 
load  transfer  was  examined  for  an  aiL  aluminum  fuseLage  panel  with  a crack  in 
the  skin  located  normal  to  the  longerons  and  parallel  to  the  frames.  Comparisons 
are  made  between  analytical  and  experimental  strain  and  residual  strength  data 
for  biaxial ly  loaded  structure, 

A summary  of  the  technique  is  presented  which  will  enable  analytical  pre- 
dictions of  amount  of  slow  tear,  crack  arrest,  and  residual  strength,  as  well 
as  defining  when  a stiffener  or  skin  critical  condition  prevails,  A compari- 
son of  the  accuracy  of  the  method  with  experimental  residual  strength  data 
shows  that  predictions  can  be  made  which  are  well  within  the  scatter  in  the 
materiai  crack  growth  resistance  data  which  is  ± 5 percent. 


UNCLASSIFIED 


SECURITY  cl  ASSt  FICATSON  OF  THIS  PAGEfWh*fi  &■*#  Bnttf+4) 


PREFACE 


This  report  was  prepared  by  Northrop  Corporation,  Aircraft  Division, 
Hawthorne,  California,  under  Air  Force  Contract  F33615-72-C-1796,  The 
project  was  initiated  under  Project  Number  486U,  Task  0204,  "Advanced 
Metallic  Structures,"  Advanced  Development  Program,  The  work  reported 
herein  was  administered  under  the  direction  of  the  Air  Force  Systems 
Command,  Wright-Patterson  Air  Force  Base,  Ohio,  by  Dr,  George  Sendeckyj 
(FBEC)  and  Captain  George  F.  Zielsdorff  (ASD),  Project  Engineers, 

The  research  was  conducted  between  January  1975  and  June  1975  as  a con- 
clusion to  Phase  III  and  the  overall  program.  This  report  was  submitted  by 
the  authors  M.  M.  Ratwani  and  D,  P,  Wilhem  in  June  1975  for  AFFDL  review. 

The  report  has  been  assigned  NOR  75-86  for  internal  control  at  Northrop 
Corporation, 

The  authors  wish  to  acknowledge  the  support  of  P,  Finwall,  J,  FitzGerald, 
and  B,  J,  Mays  for  their  assistance  during  the  structural  testing  program,  and 
to  Messrs,  V,  C,  Frost  and  P,  W,  Warren  for  assistance  in  structural  test 
specimen  design  and  construction.  The  secretarial  assistance  of  Lessie  Speciale 
is  also  acknowledged. 


iii 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I INTRODUCTION  1 

II  STIFFENED  UN LARI  ALLY  LOADED  PANELS  2 

2.1  GENERAL  INTRODUCTION  2 

2.2  BROKEN  AND  INTACT  STRINGER -THIN  SKIN  2 


2.2.1 

Description  of  Test  Panel (s) 

2 

2.2.2 

Finite  Element  Modeling  and  Analysis 

5 

2.2.2. 1 

Finite  Element  Modeling  of  Intact 
and  Broken  Stringer  Panels 

5 

2. 2.2. 2 

Elastic  and  Elastic-Plastic  Analysis 
of  Panels  with  Intact  and  Broken 

S tringer 

8 

2.2.3 

Experimental  Results  and  Discussion 

14 

2. 2.3.1 

Test  Arrangement  - General 

14 

2. 2. 3. 2 

Fatigue  Precracking,  Strain  Data  and 
Fracture  of  Broken  and  Intact  Stringer 
Cases 

18 

2. 2. 3. 3 

Comparison  of  Experimental  and 

Analytical  S trains 

28 

2.2.4 

Residual 

Strength  Prediction 

36 

2.2.4. 1 

Residual  Strength  of  Panel  wLth  Intact 
Central  Stringer 

36 

2. 2.4. 2 

Residual  Strength  of  PaneL  with 

Broken  Central  Stringer 

40 

TITANIUM  PANEL 

40 

2.3.1 

Description  of  Test  Pane i 

40 

2.3.2 

Finite  Element  Modeling  and  Analysis 

41 

2. 3. 2.1 

Finite  Element  Modeling  of  the  All 
Titanium  Panel 

43 

2.3.2. 2 

Elastic  and  Elast ic -P lastic  Analyses 

43 

of  All  Titanium  Panel 


v 


TABLE  OF  CONTENTS  (CONTINUED) 


SECTION  PAGE 

2,3,3  Experimental  Results  and  Discussion  43 

2.3.3. 1 /Sr  Crack  Growth  Resistance  Data  for  47 

Beta  Mill  Annealed  Ti-6AL-4V 

2. 3. 3. 2 Fatigue  Precracking,  Strain  Data  47 

and  Fracture  of  Titanium  Panel 

2. 3*3. 3 Comparison  of  Experimental  and  56 

Analytical  Strains 

2,3*4  Residual  Strength  Prediction  56 

2.4  INTACT  STRINGER  - THICK  SKIN  65 

2.4.1  Description  of  Test  Panel  65 

2.4.2  Finite  Element  Modeling  and  Analysis  67 

2.4.3  Experimental  Results  and  Discussion  67 

2.4,3, L Fatigue  Precracking  Strain  Data  67 

and  Fracture  of  Thick  Skin  Panel 

2. 4. 3. 2 Comparison  of  Experimental  and  78 

Analytical  Strains 

2.4.4  Residual  Strength  Prediction  78 

III  STIFFENED  B1AXIALLY  LOADED  PANELS  84 

3.1  GENERAL  INTRODUCTION  84 

3.2  DESCRIPTION  OF  TEST  PANEL(S)  35 

3.3  FINITE  ELEMENT  MODEL  AND  ANALYSIS  89 

3.3.1  Finite  Element  Modeling  of  Biaxial ly  Loaded  89 

P ane 1 s 

3.3.2  Elastic  and  Elastic -P Lastic  Analysis  of  the  89 

Biaxial  Panel 

3.3.2. 1 Elastic  Analysis  of  Biaxial  Panel  89 

3. 3. 2. 2 Elastic -P lastic  Analysis  of  99 

Biaxial  Panel 


vi 


TABLE  OF  CONTENTS  (CONTINUED) 


SECTION  PAGE 

3.4  EXPERIMENTAL  RESULTS  AND  DISCUSSION  99 

3.4.1  Fatigue  Precracking,  Strain  Data  and  Fracture  of  103 

Biaxial  Loaded  Panels 

3.4.2  Comparison  of  Experimental  and  Analytical  Data  120 

3.5  RESIDUAL  STRENGTH  PREDICTION  128 

IV  SUMMARY  OF  OVERALL  PROGRAM  AND  PREDICTION  METHODS  132 

4.1  PROGRAM  SUMMARY  132 

4.2  PREDICTIVE  PROCEDURE  (STEPS)  133 

V CONCLUSIONS  134 

5.1  GENERAL  TRENDS  134 

5.2  OVERALL  ACCURACY  OF  THE  METHOD  134 

REFERENCES  136 


vii 


LIST  OF  ILLUSTRATIONS 


F £ Ru  re  Pane 

1 Drawing  of  Wing  Panel  - Angle  Stiffened  3 

2 Basic  Finite  Element  Model  for  Angle  Stiffened*  Lower  Wing  7 

PaneL  s 

3 Contours  for  J integral  Calculations  - Angle  Stiffened*  Lower  9 

Wing  Panels 

4 Comparison  of  Square  Root  of  J VaLues  for  Intact  and  Broken  JO 

Central  Stringer  Cases  - Thin  Skin 

5 Crack  Opening  Displacement  for  Broken  (-23)  and  Intact  (-15)  II 

Stringers*  Angle  Stiffened  Panels 

6 Prediction  of  Broken  Stringer  Panel  Stress  From  Intact  12 

Stringer  Analysis 

7 Stress  in  Central  Stringer  as  a Function  of  Applied  Stress  13 

Based  on  Elastic  and  Elastic -P las  tic  Analysis  - Intact  Stringer 

8 Square  Root  of  J as  a Function  of  Applied  Stress  for  Varying  15 

Crack  Lengths  - Intact  Stringer 

9 Square  Root  of  J as  a Function  of  Crack  Length  - Intact  is 

S t ringer 

10  View  of  Test  Panel  in  500,000  Pound  Load  Frame  17 

11  Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened,  19 

Lower  (Center  Stringer  Broken) 

12  Overall  View  of  Skin  Side  of  Wing  Panel  AngLe  Stiffened,  20 

Lower  (Center  Stringer  Intact) 

13  Overall  View  of  Angle  Stiffened  Side  of  Wing  Panel  Angle  21 

Stiffened*  Lower  (Center  Stringer  Broken) 

14  Overall  View  of  Angie  Stiffened  Side  of  Wing  Panel  Angle  22 

Stiffened*  Lower  (Center  Stringer  Intact) 

15  Strain  Gage  Locations  -23  and  -15  Panels  25 

16  Load -Displacement  Curves  for  -23  Panel  Angle  Stiffened*  26 

Lower  (Center  Stringer  Broken) 


viii 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


Figu re 

Page 

17 

Load -Displacement  Curves  for  -15  Panel  Angle  Stiffened* 

Lower  (Center  Stringer  Intact) 

27 

18 

View  of  Fracture  Surface  from  Skin  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Broken) 

29 

19 

View  of  Fracture  Surface  from  Skin  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Intact) 

30 

20 

View  of  Fracture  Surface  from  Stringer  Side,  Wing  Panel 

Angle  Stiffened,  Lower  (Center  Stringer  Initially  Broken) 

31 

21 

View  of  Fracture  Surface  from  Stringer  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Intact) 

32 

22 

Strain  in  Central  Stringer  for  -15  Panel,  Angle  Stiffened, 

Lowe  r 

33 

23 

Strains  in  Central  Attached  Angle  Away  from  the  Crack  Plane 
as  a Function  of  Applied  Stress  for  -15,  Intact  Angle 

Stiffened  Panel 

34 

24 

Strain  in  Central  Stringer  for  -23  Panel,  Angle  Stiffened, 
Lower 

35 

25 

Strain  in  Central  and  Outer  Attached  Angles  as  a Function  of 
Applied  Stress  for  Broken  (-23)  and  Intact  (-15)  Stiffened 

Pane  1 s 

37 

26 

Residual  Strength  Prediction  Plot  for  -15  Panel  (Intact 
Stringer) 

39 

27 

Strain  Cage  Locations  - Titanium  Panel 

42 

28 

Elastic  and  Elastic  Plastic  Analysis  of  All  Titanium  Panel 

44 

29 

Square  Root  of  J as  a Function  of  Applied  Stress  for  Varying 
Crack  Lengths  - All  Titanium  Panel 

45 

30 

Square  Root  of  J as  a Function  of  Crack  Length  - Titanium 

Pane  l 

46 

31 

/Jg  Crack  Growth  Resistance  Data  - T1-6A1-4V  (Beta  Mill 
Annealed)  - .053M,  LT 

48 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


49 

50 

51 

52 

53 

55 

57 

58 

59 

60 

61 

62 

64 

66 

68 

69 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


/JR  Crack  Growth  Resistance  Data  - TL-6A1-4V  (Beta  Mill 
Annealed)  - .053°,  TL 

/Jr  Crack  Growth  Resistance  Data  - T1-6AL-4V  (Beta  Mill 
Annealed)  - , 1 87" f LT 

/JR  Crack  Growth  Resistance  Data  * Ti-6AL-4V  (Beta  Mill 
Annealed)  - .L87'\  TL 

Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened, 
Lower  (All  Titanium  Panel) 

Overall  View  of  Angle  Stiffened  Side  of  Wing  Panel  Angle 
Stiffened,  Lower  (All  Titanium  Panel) 

Strain  Gage  Locations  - Titanium  Panel 

Load -Di sp lac ement  Curves  for  Titanium  Wing  Panel  Angle 
Stiffened,  Lower 

View  of  Fracture  Surface  From  Skin  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (All  Titanium  Panel) 

View  of  Fracture  Surface  From  Stringer  Side,  Wing  Panel 
Angle  Stiffened,  Lower  (All  Titanium  Panel) 

Strain  in  Central  Stringer  for  All  Titanium  Panel,  Angle 
Stiffened,  Lower 

Strain  in  Central  Attached  Angle  Away  From  the  Crack  Plane 
as  a Function  of  Applied  Stress  for  All  Titanium  Panel 

Elastic  and  E las ti c -F Las t ic  Stress  in  Central  Stringer  as  a 
Function  of  Applied  Stress  - AIL  Titanium  Panel 

Residual  Strength  Prediction  Plot  for  All  Titanium  Panel 

Strain  Gage  Locations  - Thick  Skin  (-1)  Panel 

Elastic  and  E La stic-P lastic  /J  Values  for  Thick  Skin,  -I 
Panel  as  a Function  of  Applied  Stress 

Elastic  and  Elastic -P lastic  Stress  in  Central  Stringer  as  a 
Function  of  Applied  Stress  -1  Panel 


x 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


F iRure  Page 

48  Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened,  71 

Lower  (Thick  Skin) 

49  Overall  View  of  Angle  Stiffened  Side  of  Wing  Panel  Angle  72 

Stiffened,  Lower  (Thick  Skin) 

50  Load  - Displacement  Curves  for  Thick  Skin  - 7075-T73  Wing  74 

Panel  Angle  Stiffened,  Lower 

51  Load -Displacement  Curves  to  Fracture  for  Thick  Skin  - 75 

7075-T73  Wing  PaneL  Angle  Stiffened,  Lower 

52  View  of  Panel  Fracture  From  Stringer  Side,  Wing  Panel  Angle  75 

Stiffened,  Lower  (Thick  Skin) 

53  Close-up  View  of  Fractured  Angles,  Wing  Panel  Angie  77 

Stiffened,  Lower  (Thick  Skin) 

54  Strain  Gage  Locations  - Thick  Skin  (-1)  Panel  79 

55  Strain  in  Central  Stringer  for  Thick  Skin,  -1  Panel,  Angle  80 

S ti f f ened  , Lower 

56  Strain  in  Central  Stringer  as  a Function  of  Applied  Stress  81 

- 1 Panel 

57  Residual  Strength  Prediction  Plot  for  Thick  Skin,  -1  Panel  82 

58  Fuselage  Panel-Frame  and  Stringer  Stiffened  (Biaxial  Panel)  87 

59  Finite  Element  Model  of  Biaxial  Panel (s)  90 

60  Contours  Used  to  Evaluate  J for  Biaxial  Panels  91 

61  Square  Root  of  J for  Varying  Crack  Lengths  as  a Function  of  93 

Applied  Stress  - Biaxial  Panel 

62  Stress  in  the  Central  Longeron  as  a Function  of  Gross  Area  9 4 

Panel  Stress  for  Varying  Crack  Length  - Biaxial  Panel 

63  Elastic  </J  Values  as  p Function  of  Applied  Stress  for  Two  95 

Biaxial  Load  Ratios 

64  Stress  in  the  Central  Longeron  as  a Function  of  Applied  Stress  96 
for  Two  Biaxial  Load  Ratios 


xi 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


Figure 

65  Total  Crack  Surface  Displacements  as  a Function  of  Distance  97 

Along  the  Crack  for  Two  Biaxial  Load  Ratios 

66  Stresses  in  Adjacent  Longerons  for  Two  Biaxial  Load  Ratios  98 

67  Square  Root  of  J for  Varying  Biaxial  Load  Ratios  Using  100 

P rand t T - Reus s Material  Assumptions 

68  Stress  In  Central  Longeron  for  Varying  Biaxial  Load  Ratios  101 

Using  P rand 1 1 - Reuss  Material  Assumptions 

69  Plastic  Zone  Sizes  for  Three  Biaxial  Load  Ratios  Using  102 

P rand tl -Reuss  Material  Assumptions 

70  Stiffener-Side  View  of  Loading  Arrangement  for  Biaxial  109 

P anel s 

71  Side  View  of  Lateral  Loading  Arrangement  for  Biaxial  105 

Panels 

72  Skin- Side  View  of  Biaxial  Test  Panel  107 

73  Longeron  and  Frame  Side  View  of  Biaxial  Test  Panel  108 

74  Strain  Gage  Locations  - Biaxial  Panels  110 

75  Load-Di splacement  Curves  for  Biaxial  Panel  2-2  During  112 

Strain  Survey  and  Slow  Tear 

76  Load -Displacement  Curve  for  Biaxial  Panel  2-2  to  Fracture  H3 

After  Slow  Tear 

77  Load -D i splacement  Curve  for  Biaxial  Panel  2-3  to  Fracture  114 

78  Fracture  Sequence  for  2-2  Biaxia 1 1 y Loaded  Panel  (60  kips  115 

Lateral  Load) 

79  Fracture  Sequence  for  2-3  Biaxial ly  Loaded  Panel  (120  kips  118 

Lateral  Load) 

80  View  of  Fracture  Surface  From  Skin  Side,  2-2  BLaxially  121 

Loaded  Pane  1 

81  View  of  Fracture  Surface  From  Longeron /F rame  Side,  2-2  122 

Biaxially  Loaded  Panel 

82  View  of  Fracture  Surface  From  Skin  Side,  2-3  Biaxially  123 

Loaded  Panel 


xii 


LIST  OF  ILLUSTRATIONS  (CONTINUED) 


F i gu  re 

Page 

83 

View  of  Fracture  Surface  From  Longeron/F rame  Side, 

2-3  Biaxialiy  Loaded  Panel 

124 

84 

Strain  in  Central  Longeron  Away  From  the  Crack  at  Three 

Biaxial  Load  Ratios 

125 

85 

Strain  in  Longeron(s)  at  0.6  Inches  Away  from  the  Crack, 

BiaxlaL  Panel  2-2 

126 

86 

Strain  in  Frame  Member  as  a Function  of  Applied  Stress  for 
Three  Biaxial  Load  Ratios 

127 

87 

Total  Crack  Opening  Displacements  for  Two  Biaxial  Load  Ratios 
as  a Function  of  Applied  Stress 

129 

88 

Square  Root  of  J Versus  Crack  Length  Curve  for  Biaxial  Panel 

130 

i 


xiii 


LIST  OF  TABLES 


?ableB  Page 

I THIN  SKIN -ANGLE  STIFFENED  WING  PANEL  AREAS  AND  MATERIAL  6 

IDENTIFICATION 

II  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA,  -23  AND  -15  PANELS  18 

III  STRAIN  GAGE  DATA,  -23  WING  PANEL  ANGLE  STIFFENED,  LOWER  (CENTER  23 

STRINGER  BROKEN) 

IV  STRAIN  GAGE  DATA,  -15  WING  PANEL  ANGLE  STIFFENED,  LOWER  (CENTER  24 

STRINGER  INTACT) 

V TITANIUM  SKIN-ANGLE  STIFFENED  WING  PANEL  AREAS  AND  MATERIAL  41 

IDENTIFICATION 

VI  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA  - TITANIUM  PANEL  47 

VII  STRAIN  GAGE  DATA,  TITANIUM  WING  PANEL  STIFFENED,  LOWER  54 

VUI  THICK  SKIN  (-1)  - ANGLE  STIFFENED  WING  PANEL  AREAS  AND  MATERIAL  65 

IDENTIFICATION 

IX  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA  - THICK  SKIN  PANEL  70 

X STRAIN  GAGE  DATA,  -1  THICK  SKIN  WING  PANEL  ANGLE  STIFFENED,  73 

LOWER 

XI  FRAME -LONGERON  STIFFENED  PANEL  AREAS  AND  MATERIAL  1DENTIFICA-  85 

TION 

XII  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA,  2-2  AND  2-3  BIAXIAL  103 

PANELS 

XIII  STRAIN  GAGE  DATA,  2-2  BIAXIAL  PANEL  109 

XIV  STRAIN  GAGE  DATA,  2-3  BIAXIAL  PANEL  111 

XV  COMPARISON  OF  PREDICTED  AND  MEASURED  RESIDUAL  STRENGTHS  FOR  135 

PHASE  II  AND  PHASE  III  PANELS 


xiv 


I INTRODUCTION 


This  report  describes  the  research  conducted  during  Phase  ILL  of  a three- 
phase  investigation  into  the  development  of  an  improved  method  of  thin  section 
residual  strength  prediction  of  cracked  aircraft  structure  (where  conditions 
of  plane  stress  or  mixed  mode  fracture  prevail).  Phases  I and  II  summarized 
the  current  state-of-the-art  of  structural  residual  strength  prediction  and 
established  an  analytical  procedure  and  fracture  criterion  to  predict 
residual  strength  (Reference  1).  Using  a flexible  fastener  model  and  the 
Bueckner/Hayes  approach  to  Dugdale  type  elastic-plastic  analysis  good  pre- 
dictions of  residual  strength  for  both  skin  and  stringer  critical  conditions 
were  made  in  Phase  II  for  two-bay,  uniaxially  loaded  panels.  The  studies 
undertaken  during  Phase  III  were  aimed  at  defining  the  limits  of  the  technique 
as  well  as  proving  (and  improving)  its  ability  to  predict  both  crack  arrest 
and  fracture  in  multi-bay  construction  and  multi-axis  loading. 

In  Section  II  the  following  structural  test  panels  were  analyzed  and 
tested:  Thin  skin  (all  aluminum)  pane  Its)  with  both  intact  and  broken 

stringers,  a thick  skin  (all  aluminum)  panel  with  mixed  mode  fracture  behavior, 
and  an  all  titanium  panel*  Ml  loading  was  performed  in  uniaxial  tension  and 
excellent  correlation  was  obtained  with  the  experimental  data*  Both  crack 
arrest  and  stringer  critical  conditions  could  be  predicted. 

Two  biaxial ly  loaded  panels  are  discussed  in  Section  111  where  a complete 
P rand 1 1 - Reuss  analysis  was  performed  for  the  crack  length  of  test.  The 
influence  of  biaxial  stress  on'  crack  opening  displacement,  plastic  zone  size, 
and  load  transferred  to  the  stiffener  were  determined.  Prediction  of  residual 
strength  for  the  one  panel  which  failed  in  the  test  area  were  made  based  on 
certain  approximations  with  a good  degree  of  success.  The  unknowns  associated 
with  the  biaxial  failure  criterion  were  also  examined  and  means  proposed  to 
study  load  interaction  effects, 

A summary  of  the  overall  three-phase  program  effort  is  included  in 
Section  IV,  which  contains  the  ten  steps  required  to  predict  structural 
residual  strength. 

In  Section  V,  a short  summary  is  presented  of  the  conclusions  to  the 
Phase  HI  research.  Also  included  is  a table  which  compares  analytical  and 
experimental  results  for  both  phases  of  the  program  which  indicates  an 
average  difference  of  approximately  four  percent  in  predicted  versus  measured 
residual  strength  for  all  panels  tested, 

it  must  again  be  noted  that  the  proposed  method  of  analysis  given  in 
Section  4,2  is  not  intended  for  use  in  performing  parametric  studies  (i.e., 
stiffener  spacing,  rivet  pitch,  etc,)  which  could  be  more  cheaply  accomplished 
using  linear  elastic  assumption  and  closed  form  solutions*  However,  the 
proposed  method  does  offer  the  ability  to  predict  residual  strength  accurately 
for  a “fatigue  sized"  structure  which  is  expected  to  fracture  with  large 
amounts  of  slow,  stable  tear  and  plasticity  prior  to  structural  failure. 
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II  STIFFENED  UNI  AXIALLY  LOADED  PANELS 


To  demonstrate  the  utility  of  the  residual  strength  prediction  method 
for  those  structures  loaded  primarily  in  tension  and  for  different  materials 
and  skin  thickness  four  panels  were  analyzed,  fabricated,  and  tested.  This 
Section  describes  those  analyses,  panels,  and  test  data. 

2.1  GENERAL  INTRODUCTION 

The  residual  strength  predictive  method  should  be  able  to  predict  fracture 
in  panels  of  typical  aircraft  structural  geometry  and  materials*  For  this 
reason  it  was  elected  to  examine  lower  wing  structure  since  high,  uniaxial 
stress  occurs  in  this  region*  Skin  material  was  7075-173  of  two  thicknesses, 
0,063  and  0,193  inch,  selected  to  produce  essentially  plane  stress  and  mixed 
mode  fracture  behavior.  In  addition,  a panel  consisting  of  all  titanium 
material , sized  to  represent  typical  fighter  aircraft  lower  wing  structure, 
was  analyzed  and  tested  using  the  proposed  residual  strength  prediction 
technique.  The  reinforcing  substructure  for  all  panels  consisted  of  angle 
sections  which  were  fastened  to  the  skin  using  standard,  aircraft  fastening 
techniques.  For  one  panel  configuration  both  an  initially  broken  and  intact 
central  stiffener  were  analyzed  and  tested. 

The  analytical  techniques  employed  in  predicting  the  residual  strength 
of  the  uniaxial  loaded  panels  of  this  Section  were  those  outlined  in  Reference 
1 which  consists  of  employing  the  Bueckner/Hayes  energy  approach  using 
Dugdale  plasticity  assumptions  combined  with  a flexible  fastener  finite 
element  model.  Failure  is  based  on  either  a stiffener  or  skin  critical 
failure  criterion.  The  basis  of  the  skin  critical  criterion  is  the  tangency 
between  the  material  crack  growth  resistance  curve,  or /Jg  curve,  with  the 
crack  driving  force  curves  in  terms  of  /J, 

2.2  BROKEN  AND  INTACT  STRINGER -THIN  SKIN 


The  panel  geometry  selected  for  this  analysis  consists  of  typical  lower 
wing  construction  (outboard  of  the  engines)  for  a large  transport  type  aircraft* 
Crack  geometries  were  similar  for  two  identical  paneLs.  The  difference  between 
the  two  was,  one  contained  an  initially  broken  central  stringer  and  it  was 
intact  in  the  other.  However,  It  was  realized  that  either  this  situation  or 
the  through  cracked  skin- stringer  could  occur  in  service  situations  in  a 
random  manner.  Therefore,  the  panels  analyzed  and  tested  In  this  part  of 
Phase  III  represent  those  situations  where  a common  crack  length  was  used  for 
comparative  purposes  in  a broken  and  intact  stringer  situation, 

2,2*1  Description  of  Test  FaneKs) 

The  drawing  of  the  test  panel  (s)  is  shown  in  Figure  1 and  is  identified 
as  -15  assembly  for  the  intact  stringer  and  -23  for  the  broken  central  stringer 
case*  The  overall  skin  width  of  the  7075-T73  was  38*5  inches  with  a 5.5  inch 
spacing  between  rivet  and  angle  centerLines,  This  spacing  produced  a six  bay 
paneL,  symmetric  about  the  panel  centerline,  i.e*,  three  bays  on  either  side 
of  the  center  angLe* 
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The  selected  angle  sections  were  AND  10134*1001  (7073-T6  extrusions) 
with  the  one  Inch  leg  attached  by  3/16  inch  diameter*  flush  head  steel  HI-LOK's 
(19PB-3)  at  1.3  Inch  pitch  in  the  test  area.  In  that  portion  of  the  lower 
wing  simulated  by  these  panels*  the  upstanding  leg  would  normally  be  attached 
to  the  upper  wing  through  an  attached  shear  web  (chordwise).  Table  I lists 
the  panel*  skin  and  angle  cross  sectional  thicknesses*  areas*  and  materials 
(including  skin  material  identification,  see  Reference  2 for  code).  It  should 
be  noted  that  the  29%  stiffening  ratio  in  Table  I is  representative  of  con* 
stmction  for  typical*  large  transport  outer  lower  wings. 

Figure  1 shows  that  the  initial  starter  slot  had  a half  length  of  2.50 
inches  and  an  overall  length  of  3,00  inches  for  both  the  -L5  and  -23  panels. 
Overall  panel  lengths  were  normally  96  inches  including  padded  grip  ends  and 
74  inches  excluding  grip  ends.  The  padded  ends  were  designed  to  produce  a 
minimal  amount  of  inplane  bending  due  to  loading  eccentricity.  Each  panel 
contained  20  strain  gages*  the  location  of  which  are  detailed  in  Section 
2.2,3,  A standard  beam  type  clip  gage  was  mounted  on  the  skin  side  of  both 
panels  to  record  crack  opening  displacemnt  at  the  crack  and  panel  centerline. 

2,2,2  Finite  Element  Modeling  and  Analysis 

The  analyses  of  the  panels  with  intact  (-15)  and  broken  (-23)  central 
angle  were  conducted  in  a manner  similar  to  the  zee  stiffened  panels  of 
Phase  II  (Reference  1)  of  this  study.  Both  elastic  and  elastic-plastic  analyse 
were  conducted  for  various  crack  lengths  for  the  panel  with  the  intact  stringer 
The  elastic  - plastic  analysis  was  based  on  the  assumption  of  a Dugdale  type 
strip  plastic  zone.  For  the  panel  with  the  broken  stringer*  elastic  and 
elastic-plastic  analyses  were  conducted  for  the  crack  length  used  in  the 
test  panel.  The  finite  element  modeling*  elastic,  and  elastic  - plastic 
analyses  of  the  panels  are  discussed  in  the  following  subsections, 

2 , 2 , 2 , 1 Finite  Element  Modeling  of  Intact  and  Broken  Stringer  Panels 

The  finite  element  modeling  is  similar  to  that  used  in  Phase  IL  (see 
Reference  1).  The  panel  is  modeled  as  a two-dimensional  structure.  Triangu- 
lar membrane  elements  were  used  to  model  the  skin.  The  angle  stiffeners  of 
this  panel  are  modeled  Like  the  zee  stiffeners  for  the  zee  stiffened  panels 
of  Phase  II.  The  connected  leg  of  the  angle  stiffener  is  modeled  as  a rod 
element  and  the  upstanding  leg  of  the  angle  is  modeled  as  a rectangular 
membrane  element.  The  effect  of  asymmetry  caused  by  the  L shaped  cross- 
section  is  not  considered.  In  this  case  only  a quarter  panel  need  be  modeled 
for  finite  element  analysis  due  to  panel  symmetry-  The  analysis*  of  Phase  II 
panels  showed  that  in  cracked*  stiffened  panels  only  four  to  eight  fasteners 
above  (and  below)  the  crack  surface  connecting  stringer  to  skin  were  effective 
in  load  shedding  from  skin  to  stringer.  The  other  fasteners  had  little 
Influence  on  load  transfer  to  the  stringers.  The  number  of  fasteners  effective 
in  Load  transfer  was  also  found  to  be  dependent  on  crack  length.  Taking  into 
consideration  these  trends  only  the  first  nine  fasteners  were  modeled  as 
individual  shear  eLements  in  the  finite  element  model.  The  depth  of  the  shear 
elements  was  kept  the  same  as  the  diameter  of  the  fastener.  That  portion  of 
the  panel  beyond  the  nine  fasteners  from  crack  plane  was  assumed  to  be  con- 
tinuously attached.  This  assumption  allows  a coarser  grid  to  be  used  in  the 


5 


TABLE  I THIN  SKIN -ANGLE  STIFFENED  WING  PANEL 
AREAS  AND  MATERIAL  IDENTIFICATION 


PANEL 

NUMBER 

MEAS 

THLC 

(In 

URED 

KNESS 

c he  s ) 

CROSS 

SECTIONAL 

AREAS 

( Inches^) 

SKIN 

MATERIAL/ 

AND  l.D. 

(See  Ref.  2) 

PERCENTAGE 

OF  TOTAL 
ANGLE  AREA 
TO  SKIN 

AREA 

Angle 

Skin 

Total 

Aneles 

Skin 

-15 

.063 

.0625 

.700 

2.40625 

7075-T  73/4-3 

29.17= 

TOTAL  PANEL 

CROSS 

SECTIONAL 

AREA  = 3,10625 

-23* 

.063 

.0625 

,700 

2.40625 

7075-T73/4-2 

29.  17. 

TOTAL  PANEL 

CROSS 

SECTIONAL 

AREA  = 3.10625 

^Central  Stringer  Broken 


finite  element  modeling  away  from  the  crack , reducing  the  total  number  of  grid 
points  and  computer  run  time.  This  is  a particularly  important  consideration 
in  this  type  of  panel  where  four  stiffeners  are  modeled  in  one  quarter  panel. 
If  all  the  fasteners  in  alL  the  stiffeners  are  modeled  as  individual  shear 
elements,  the  number  of  grid  points  becomes  excessive  with  Large  computer 
run  times.  The  flexibLe  fastener  model  was  used  in  the  analysis.  The  shear 
elements  in  the  finite  element  model  were  proportioned  to  account  for  the 
flexibility  of  the  fasteners  as  discussed  in  Section  3.1.3  of  Reference  1, 

The  finite  element  model  of  the  angle  stiffened  panel  is  shown  in  Figure  2* 

The  model  has  710  grid  points  and  contains  781  triangular  membrane  elements 
representing  the  skin.  The  size  of  the  smallest  triangular  element  used 
ahead  of  the  crack  tip  is  O.L25  inch.  The  computer  run  time  for  an  elastic 
analysis  with  a crack  length  of  5.50  inches  is  5 minutes,  A compLete  Dugdale 
model  type  elastic- plastic  analysis  for  the  same  crack  length  with  eight  load 
(or  plastic  zone)  cases  requires  approximately  22  minutes. 
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Figure  2.  Basic  Finite  Element  Model  for  Angle  Stiffened*  Lower  Wing  Panel 


2 . 2 * 2 , 2 Elastic  and  E lastic -P last ic  Analysis  of  Panels  With  Intact  and 
Broken  Stringer 


An  elastic  analysis  of  the  intact  1-15)  and  broken  (-23)  central  stringer 
panels  was  conducted  for  a half  crack  length  of  2,75  inches.  Two  different 
contours  were  used  to  compute  J integral  values.  These  two  contours  are 
shown  in  Figure  3,  The  elastic  square  root  of  J (/J)  values  obtained  for  the 
two  contours  were  within  0,5  percent.  Figure  4 shows  the  plot  of  square 
root  of  J versus  appLied  stress  normalized  to  skin  material  yield  strength, 
for  these  two  panels.  It  is  seen  that  the  square  root  of  J values  for  the 
broken  central  stringer  case  is  approximately  22  percent  higher  than  those 
for  the  intact  stringer  case.  The  crack  surface  openings  for  these  two  cases 
are  plotted  in  Figure  5.  The  crack  surface  opening  at  the  center  line  of  the 
crack  in  the  broken  stringer  case  is  approximately  77  percent  higher  than 
that  for  the  intact  central  stringer  panel.  It  should  be  noted  in  Figure  5, 
that  the  maximum  crack  surface  opening  for  the  intact  stringer  case  occurs 
away  from  the  center  line  of  the  crack,  whereas  for  the  broken  stringer  case 
it  occurs  at  the  center  line  of  the  crack. 

In  the  elastic-plastic  analysis  of  the  -15  and  -23  panels  a Dugdale  type 
strip  plastic  zone  was  assumed  ahead  of  the  crack  tip*  The  Bueckne r-Ha yes 
approach,  as  discussed  in  Reference  1 , was  used  for  the  Dugdale  type  elastic- 
plastic  analysis.  The  variation  of  /J  with  p/Ft^  for  both  intact  (-15)  and 

broken  (-23)  central  stringer  cases  using  this  analysis  is  shown  in  Figure  6 
for  a half  crack  length  of  2,75  inches.  For  the  same  applied  stress  the 
square  root  J values  for  the  broken  stringer  case  are  higher  than  those  for 
the  intact  stringer  case.  At  an  applied  stress  of  80  percent  of  the  skin 
material  yield  the  broken  stringer  case  /J  value  is  22  percent  higher  than 
that  for  an  Intact  stringer  case*  For  both  panels  the  elastic -plastic  analysis 
gives  square  root  J values  higher  than  those  given  by  an  elastic  analysis. 

The  difference  in  the  /J  values  between  the  two  methods  of  %/J  analyses  is 
12  percent  for  an  applied  stress  of  80  percent  of  the  skin  material  yield 
strength.  For  the  intact  central  stringer  case  the  variation  of  stress  in  the 
central  stringer,  at  the  centerline  of  the  crack  and  panel  with  applied  stress 
is  shown  in  Figure  7.  At  an  applied  stress  of  80  percent  of  yield  the  stresses 
in  the  stringer  from  the  e lastic - plas t ic  analysis  are  20  percent  higher  than 
those  given  by  a purely  elastic  analysis.  These  stresses  have  been  plotted  to 
determine  if  the  crack  length  under  consideration  represents  a stringer 
critical  case  as  proposed  in  the  analysis  and  discussed  subsequently  in 
Sub-Section  2,2.4. 

In  the  residual  strength  prediction  method  it  is  necessary  to  plot  a 
curve  between  the  square  root  J and  crack  length  for  various  applied  loads. 

In  order  to  determine  these  values  it  is  necessary  to  perform  an  elastic- 
plastic  analysis  for  various  crack  lengths.  The  selected  half  crack  lengths 
were  1,50,  2,75,  4.0,  5.0,  5.5,  and  6,5  inches  for  the  intact  stringer  case 
(-15  panel)  only,  A similar  analysis  was  not  conducted  for  the  broken 
stringer  (-23  panel)  case  due  to  the  large  computer  run  times  involved. 

The  residual  strength  predicted  for  the  intact  stringer  case  with  the 
square  root  of  J curve  shown  in  Figure  4 were  used  to  predict  the  residual 
strength  of  the  broken  stringer  case  as  discussed  in  Sub-Section  2.2.4.  The 
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Figure  6.  Prediction  of  Broken  Stringer  PaneL  Stress  From  Intact  Stringer  AnaLysls 
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Figure  7.  Stress  In  Central  Stringer  as  a Function  of  Applied  Stress  Based 
on  Elastic  and  Elastic-Plastic  Analysis  - Intact  Stringer 
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square  root  J values  as  a function  of  applied  stress  for  various  crack  lengths 
are  shown  in  Figure  8,  These  values  are  cross  plotted  to  show  the 
variation  of  square  root  of  J with  crack  length  for  the  same  applied  stress 
in  Figure  9 * For  the  same  applied  stress,  square  root  J values  increase 
with  crack  length  to  a half  crack  length  of  about  4,0  inches.  With  further 
increase  in  crack  length  the  square  root  J values  decrease  due  to  the 
influence  of  the  adjacent  stringer.  Once  the  crack  is  beyond  the  center  line 
of  the  stringer  the  value  of  square  root  J increases  again  with  crack  length. 

2.2.3  Experimental  Results  and  Discussion 

To  prevent  repetition  within  the  various  sections,  a complete  description 
of  overall  test  arrangement,  data  recording  instrumentation  and  reduction 
technique  are  described  in  this  sub~section.  Following  this  discussion 
the  experimental  results  for  the  intact  and  broken  stringer  cases  will  be 
presented,  along  with  discussion  of  the  results  obtained. 

2. 2. 3 * I Test  Arrangement  - General 


All  panels  were  loaded  uniaxial iy  in  a 500,000  pound  load  frame.  The 
overall  test  setup  is  shown  in  Figure  10.  Loading  of  the  test  panel  is 
accomplished  through  a slotted  clevis  attached  directly  to  the  bolted  on  grip 
ends . Load  is  introduced  through  an  M , T ,S  , Mode  1 406  Controller  which  controls 
the  500,000  pound  cylinder.  Load  was  recorded  on  three  separate  units i a 
HONEYWELL  194  chart  recorder,  HEWLETT-PACKARD,  model  704 5A , x-y  recorder  and 
on  magnetic  tape  in  conjunction  with  strain  recording. 

High  speed,  16  mm  (split  frame)  black  and  white  motion  picture  coverage 
was  used  to  monitor  both  slow  stable  tear  and  final  fracture  at  approximately 
1,000  frames  per  second.  Load  was  monitored  on  film  by  placing  a DC  voltmeter 
reading  output  load  voltage,  wi thin  the  field  of  view  of  the  camera. 

Instrumentation  on  all  wing  panels  consisted  of  an  M.T.S.  clip  gage 
excited  by  an  ECTR0N  Differential  DC  Amplifier  (Model  687)#  The  output  was 
recorded  on  an  x-y  plotter.  This  gage  in  all  cases  spanned  the  centerline  of 
the  crack  and  panel  on  the  skin  side.  Knife  edges  were  mounted  on  the  panel 
with  epoxy  and  the  clip  gage  positioned  between  the  two  knife  edges. 

Recording  of  the  twenty  (20)  strain  gages  on  each  panel  was  accomplished 
in  two  ways:  first  by  direct  reading  alL  twenty  gages  employing  a DAT RAN 
Model  strain  gage  printer  during  load  holding  procedures  and  second  through 
magnetic  tape  recording  of  L2  selected  gages  during  both  load  holding  and 
increasing  Load  testing.  Three  (3)  gages  were  alternately  monitored  using 
both  methods  of  recording  to  determine  any  error  in  the  strain  sensing 
systems . 

The  tape  data  recording  system  consisted  of  the  strain  gage  signal  trans- 
mitted through  a balance  box  fed  to  an  ASTRO -DATA  Model  889  Amplifier  and 
recorded  on  a 14  channel  HONEYWELL  Model  7600  recorder.  Twelve  channels  were 
utilized  for  strain  gage  recording,  one  for  load  cell  voltage  and  the  remaining 
channel  for  voice.  After  each  test  these  data  were  played  back  through  a BRUSH 
Mark  200  recorder  for  data  reduction  of  strain  and  load  versus  time. 
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Figure  9.  Square  Root  of  J as  a Function  of  Crack  Length  - Intact  Stringer 
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leftt  Setup 


Figure  10.  View  of  Test  Panel  in  500*000  Pound  Load  Frame 
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Balancing  of  strain  gages  and  zero  initial  readings  were  taken  with  the 
specimen  hanging  "free11  from  the  top  clevis  - the  bottom  pin  was  "palled" 
during  all  initial  zero  readings*  In- test  and  post  fracture,  zero  load  strain 
readings  were  taken  with  the  Loading  pin  in  place. 

2 . 2 . 3 . 2 Fatigue  Precracking,  Strain  Data  and  Fracture  of  Broken  and 
Intact  Stringer  Cases 

In  order  to  start  from  a natural  crack  condition  in  the  two  panels,  each 
was  fatigue  precracked  approximately  0,2  inches  from  the  ends  of  the  machined 
slot  using  tension- tension , constant  amplitude,  sinusoidal  loading.  Maximum 
fatigue  cycling  load  (stress)  levels  were  determined  using  the  stress  intensity 
solutions  of  Reference  3 and  assuming  a maximum  stress  intensity  of  15  ksi 
/ inch  in  both  cases,  and  a fatigue  stress  ratio  of  approximately  0.1,  Table 
II  summarizes  the  fatigue  loadings  and  total  crack  lengths  for  these  two 
panels.  The  crack  Lengths  were  measured  after  panel  failure. 


TABLE  XI  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA s -23  AND  -15  PANELS 


FATIGUE 

LOADS 

FATIGUE  CRACK  LENGTH (S) 

MAXIMUM 

MINIMUM 

( Inches) 

PANEL 

(kips) 

(kips) 

NO.  OF  CYCLES 

LEFT 

RIGHT 

TOTAL 

-23* 

7.3 

0.8 

63,000 

2.70 

2.76 

5.46 

-13 

23.0 

3.0 

10,200 

2.79 

2.76 

5.55 

* Broken  Central  Stringer 


Figures  11  and  12  show  the  front  (skin)  side  view  of  the  -23  and  -15 
panels  in  the  Loading  frame  and  Figures  13  and  14,  the  back  (angle)  side  of 
these  paneLs,  The  positioning  of  displacement  (clip)  gage,  strain  gages 
initial  saw  slots  and  panel  detail  can  be  noted  for  each  case. 

After  fatigue  cycling  and  zero  balancing  of  all  strain  gages  and  dis- 
placement gage,  a strain/ load  survey  was  accomplished  on  each  panel  by  Load- 
ing in  specific  load  intervals.  These  data  are  noted  in  TabLe  I I I for  the 
broken  central  stiffened  panel  and  Table  IV  for  the  intact  panel.  The 
reference  of  each  strain  gage  is  shown  in  Figure  15. 

Load  versus  total  displacement  gage  output  is  shown  in  Figure  L6  for 
the  initially  broken  central  angle  panel  (-23)  and  in  Figure  17  for  the 
initially  intact  stringer  case  (-15).  For  the  -23  panel  inplane  buckling  in 
the  two  bays  was  quite  evident  at  very  low  load  levels,  (approximately  15  to 
20  kips);  whereas,  the  intact  stringer  case  required  approximately  twice  this 
load  for  crack  buckling. 
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Figure  1 L - Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened, 
Lower  (Center  Stringer  Broken) 
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Figure  L2,  Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened, 
Lower  (Center  Stringer  Intact) 
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TABLE  III  STRAIN  GAGE  DATA,  -23  WING  PANEL  ANGLE  STIFFENED,  LOWER  (CENTER  STRINGER  BROKEN) 
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TABLE  IV  STRAIN  GAGE  DATA,  -15  WING  PANEL  ANGLE  STIFFENED,  LOWER  (CENTER  STRINGER  INTACT) 
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Total  Displacement,  <3  (inches) 
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Total  Displacement , 6 (inches) 

Figure  17.  Load -Displacement  Curves  for  -15  Panel  Angle  Stiffened,  Lower  (Center  Stringer  Intact) 


Both  panels  were  loaded  in  increments  as  noted  in  the  data  of  Tables  ILL 
and  IV  and  Figures  16  and  17-  Loading  was  stopped  at  the  first  sign  of  crack 
extension  accompanying  large  plastic  zone  build-up.  After  unloading  and 
measurement  of  the  crack,  the  loading  was  re- initiated  using  ramp  loading  of 
100,000  pounds  per  minute.  Slow,  stable  tear  occurred  for  both  panels  and 
the  cracks  arrested  at  the  adjacent  stringers*  Strain  readings  were  taken  at 
the  maximum  load  point  where  the  crack(s)  arrested.  This  load  was  91  kips 

for  the  -23  panel  and  111-7  kips  for  the  -15  panel.  At  this  point  the  panels 

were  unloaded  and  a ramp  loading  of  100,000  pounds  per  minute  was  again 
initiated.  Several  load  hold  conditions  were  instituted  on  the  way  to 
fracture  at  which  time  strain  data  were  recorded.  Continuing  loading  caused 

fracture  at  107.2  kips  for  the  -23  panel  and  at  116,1  kips  for  the  -15  panel. 

Motion  picture  coverage  was  initiated  at  the  start  of  slow  tear  and  as 
close  to  fracture  as  possible  in  each  case. 

The  displacement  data  for  the  -23  panel  (Figure  16)  indicates  a stair- 
case appearance  due  to  deliberate  offsetting  of  the  x axis  zero  during  the 
initial  strain  survey  at  each  load  hold  condition.  This  offsetting  procedure 
was  not  followed  for  the  -15  panel,  except  for  the  final  fracture  trace.  On 
this  panel  (Figure  17)  the  gage  was  removed  from  the  specimen  prior  to 
fracture  since  it  was  close  to  being  out  of  travel  at  the  final  load  hold 
condition. 

Figures  18  and  19  indicate  the  fracture  surfaces  for  the  -23  and  -15 
panels.  Fracture  occurred  through  the  center  line  of  fasteners  in  both 
cases.  The  back  side  or  angle  side  of  the  -23  and  -15  panel  after  fracture 
are  shown  in  Figures  20  and  21.  The  angles  fractured  at  the  centerline 
fasteners  along  the  crack  path  in  both  cases, 

2 , 2 , 3 , 3 Comparison  of  Experimental  and  Analytical  Strains 

Figure  15  shows  the  location  of  strain  gages  on  both  the  intact  and 
broken  central  angle  panel.  Tabulated  in  Table  LV  are  the  strain  gage 
readings  for  various  applied  loads  at  a crack  length  of  5,55  inches.  Figure 
22  shows  the  experimental  and  analytical  variation  of  strain  in  the  cenLral 
stringer  with  increasing  distance  from  the  crack  for  the  -15,  intact  stringer 
panel  with  a half  c rack  length  of  2,775  inches.  Good  correlation  is  noted 
between  experimental  strains  and  those  obtained  from  the  elastic  analysis. 

The  experimental  data  of  Figure  22  is  cross  plotted  in  Figure  23  to  show 
the  variation  of  strain  with  applied  stress  for  the  -15  panel.  Note  that 
up  to  an  applied  stress  of  22  ksi  a good  correlation  is  obtained  between 
elastic  analysis  and  experimental  strains.  At  higher  applied  stresses  the 
experimental  strain  in  the  stringer  close  to  the  crack  plane  is  slightly 
Lower  than  that  given  by  an  elastic  analysis. 

For  the  broken  central  stringer  tangle)  panel  the  variation  of  strains 
in  the  centrai  stringer  with  increasing  distance  from  crack  plane  is  shown  in 
Figure  24.  It  can  be  seen  that  the  strain  gage  just  above  the  crack  surface 
(one -half  inch  away)  is  in  compression  under  tensile  loads  as  predicted  by 
the  analysis.  Good  correlation  is  obtained  between  the  elastic  analyses  and 
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Figure  18.  View  of  Fracture  Surface  from  Skin  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Broken) 
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View  of  Fracture  Surface  from  Skin  Side*  Wing  PaneL  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Intact) 
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Figure  20.  View  of  Fracture  Surface  from  Stringer  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Center  Stringer  Initially  Broken) 
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Figure  21*  View  of  Fracture  Surface  from  Stringer  Side,  Wing  Panel  Angie 
Stiffened,  Lower  (Center  Stringer  Initially  Intact) 


Figure  22,  Strain  In  Central  Stringer  for  -1.5  Panel,  Angle  Stiffened,  Lower 
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Figure  23.  Strains  in  Central  Attached  Angle  Away  from  the  Crack  Plane  as  a 
Function  of  Applied  Stress  for  -15,  Intact  Angle  Stiffened  Panel 
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Strain  In  Central  Stringer,  (Inches/Inch  x 10^) 
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the  experimental  measured  strains  up  to  an  applied  load  of  50  kips  (applied 
stress  of  16.6  ksi},  however,  at  higher  loads  the  elastic  strains  are  con- 
siderably larger.  In  comparing  Figures  22  and  24  } it  can  be  observed  that 
the  experimental  strains  deviate  from  the  elastic  results  at  a lower  applied 
stress  for  the  broken  stringer  case.  This  is  due  to  the  higher  stress 
intensity  factors  and  larger  plastic  zones  for  the  broken  stringer  panel  (-23) 
compared  to  the  intact  stringer  (-15)  panel,  A comparison  of  stringer 
strains  in  the  -15  and  -23  panels  is  shown  in  Figure  25  , The  stringer 
strains  are  from  a location  one-half  inch  above  the  crack  surface.  For  the 
broken  central  stringer  panel  (-23)  the  strain  in  the  center  stringer  is 
compressive  for  all  applied  loads  (also  see  Figure  24).  The  stringers 
located  5.5  inches  from  the  center  line  of  the  crack  and  panel  have  higher 
strains  for  the  broken  stringer  panel  compared  to  corresponding  strains  in 
the  adjacent  stringers  of  the  -15  panel.  The  correlation  between  experimental 
strains  and  elastic  analytical  strains  is  good  up  to  an  applied  stress  of 
20  ksi  for  both  panels.  Beyond  that  point  plastic  zone  bulld-up  and  small 
crack  extension  has  occurred  for  the  -15  panel. 

2.2.4  Residual  Strength  Prediction 

Ihe  residual  strength  of  panels  with  intact  and  broken  central  stringers 
ts  predicted  in  this  section.  The  results  are  compared  with  those  obtained 
from  the  tested  panels. 

2 , 2 , 4 . 1 Residual  Strength  of  Fanel  with  Intact  Central  Stringer 

The  procedure  discussed  in  the  Phase  II  report  (Reference  1)  was  used  to 
predict  the  residual  strength  of  the  complex  aircraft  structural  panels  of 
Phase  III.  The  elastic-plastic  analysis  based  on  an  assumed  Dugdale  type 
plastic  zone  ahead  of  the  crack  tip  in  conjunction  with  the  /T^  resistance 

curve  of  the  skin  material  was  used  in  the  residual  strength  prediction  for  the 
skin  critical  fracture  cases.  On  the  /J  versus  crack  size  (a)  plots  for  the 
panel,  discussed  in  Sub-Section  2.2.2, 2 (Figure  9),  the  /X  versus  aaDUV 

material  resistance  curve  is  overlayed  at  the  initial  crack  length  of  interest. 
The  applied  stress  at  which  versus  Aap^y  # and  / J versus  ^PHYSICAL  curves 

become  tangential  will  determine  the  stress  at  which  instability  occurs  and 
the  crack  propagates  rapidly.  After  this  point  of  instability,  if  the  /J 
curve  for  the  panel  dips  below  the  resistance  curve  of  the  material  the 

crack  will  be  arrested  otherwise  it  will  propagate  catastrophically.  The 
stresses  in  the  stringer  which  are  obtained  from  the  elastic  - plast ic  analysis 
and  the  ultimate  strength  of  the  stringer  material  are  used  initially  to 
determine  if  a stringer  critical  case  prevails, 

ihe  stresses  in  the  central  stringer  for  the  panel  with  the  intact  central 
stringer  and  a half  crack  length  of  2,75  inches  are  shown  in  Figure  7. 

Ihe  ultimate  strength  of  the  stringer  material  (7075-T6  extrusions)  was  found 
to  be  87.2  ksi  from  tensile  test  data.  Using  this  ultimate  stress  value 
and  the  Dugdale  type  e las  tic  - plastic  analysis  shown  in  Figure  7,  the  failure 
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of  the  central  stringer  tangle)  is  predicted  at  p/Ffc  = 0.7,  Next  consider  the 
A versus  crack  Length  plot  for  various  applied  stresses  as  shown  in  Figure  26. 
Superimposed  on  this  plot  is  the  /Jr  versus  Aap^y  resistance  curve  of  the 

7075-T73  skin  material  plotted  at  a physicaL  half  crack  Length  of  2.775  inches. 
It  is  seen  from  this  Figure  that  at  the  stringer  critical  stress  = 0,7, 

there  will  be  considerable  slow  tear  in  the  skin  and  therefore  the  stringer 
stresses  shown  in  Figure  7 wilL  not  be  valid.  From  Figure  26,  the  resistance 
curve  of  the  material  pLotted  on  the  physical  half  crack  Length  of  2.7  75 
inches  (initial  half  crack  length  in  the  tested  panel)  is  tangential  to  the 
A versus  a curve  of  the  paneL  at  an  applied  stress  of  p/F  ^ = 0,545,  Thus 

the  first  point  of  crack  instability  occurs  (after  slow  crack  growth)  at  this 
stress  and  the  crack  starts  a rapid  advance.  However,  from  Figure  26,  it  is 
noted  that  at  a half  crack  Length  of  4.45  inches  the  j~Z  curve  for  the  panel 
drops  lower  than  the  resistance  curve  of  the  material,  l.e.,  beyond  this 
point  the  resistance  of  the  material  is  higher  than  /J  developed  in  the 
stressed  panel,  and  hence  the  running  crack  will  become  arrested.  The  crack 
was  arrested  at  the  rivet  hole  in  the  -15  panel  where  the  stringer  is 
connected  to  the  skin.  The  resistance  curve  of  the  skin  material  is  now 
replotted  at  a half  crack  length  of  5.5  inches  (distance  from  centerline  of 
the  stringer  to  panel  centerline),  where  the  crack  became  arrested.  For  this 
crack  Length  the  /J-  curve  of  the  panel  becomes  tangential  to  the  resistance 
curve  at  an  applied  stress  of  p/Fty  - 0.645,  Hence  the  crack  instability 
occurs  at  that  stress  and  crack  starts  running  catast rophical Ly,  Beyond  this 
point  of  instability  the  /J  values  in  the  panel  are  higher  than  the  resistance 
curve  of  the  skin  material  and  therefore  no  possibility  exists  for  crack 
arrest.  The  /J  curve  of  the  panel  wilL  continue  to  rise  under  increasing 
load  until  the  influence  of  the  next  stringer  on  panel  stress  (11  inches  from 
the  centerline  of  the  panels  is  felt  and  the  / J value  will  once  again  have  a 
decreasing  trend  and  reach  a second  minima  at  that  point.  From  the  trend  of 
the  /J  curve  in  Figure  26  it  is  evident  that  at  an  applied  stress  of  p/F^y  = 

0,645  the  fj  value  of  the  panel  will  be  higher  beyond  a half  crack  Length  of 
7,5  inches.  Therefore  no  possibility  of  crack  arrest  is  possibLe  at  the 
second  stringer,  i.e,,  11  inches  from  the  centerline  of  the  crack,  under 
increasing  load  conditions. 

The  panel  analyzed  above  (-15  panel)  was  tested  to  failure.  Slow  tear 
in  this  panel  started  at  a load  of  60  kips  which  corresponds  to  p/Fj-y  of  0,32, 

At  this  applied  load  slow  tear  on  each  crack  tip  was  about  0,02  inch  (see 
Table  IV),  The  first  point  of  instability  (rapid  crack  extension)  occurred 
after  an  applied  load  of  60  kips  and  before  a load  of  111.7  kips.  At  111.7 
kips  the  crack  had  already  reached  the  angle  stiffeners  and  became  arrested. 

From  theory  the  first  point  of  instability  is  at  an  applied  stress  of  0,545 
F t y or  a load  of  104  kips  (see  Figure  26),  After  the  crack  had  torn  to  the 

angle  stiffeners,  the  panel  was  able  to  take  a load  of  111,7  kips  without 
any  further  slow  tear.  At  an  applied  load  of  111.7  kips  (L,e.,  p/F Ly  - 0.58) 

the  slow  tear  analysis  from  Figure  26  is  approximately  0,05  inches.  This  is 
less  than  the  radius  of  the  rivet  hole. 
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jre  26.  Residual  Strength  Prediction  Plot  for  -15  Panel  (Intact  Stringer) 
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The  measured  fai-Lure  l.oad  of  the  panel  was  116. 1 kips  (see  TabLe  IV). 

The  predicted  failure  is  at  p/Ft  = 0.645  or  panel  load  of  123.4  kips. 

This  predicted  failure  load  is  within  6.28  percent  of  the  measured  failure 
load  for  the  -15  panel. 

2. 2. 4. 2 Residual  Strength  of  Panel  with  Broken  Central  Stringer 

For  this  panel  the  central  stringer  (angle)  is  broken  and  therefore 
represents  a skin  critical  case.  A complete  elastic-plastic  analysis  for 
various  crack  lengths  was  not  conducted  for  this  panel.  An  elastic-plastic 
analysis  was  computed  for  the  tested  crack  length  of  approximately  5.5  inches. 
The  comparison  of  elastic-plastic  analysis  for  the  intact  and  broken  stringer 
panels  is  shown  in  Figure  6.  The  elastic  and  elastic-plastic  ft  values 
for  the  broken  central  stringer  case  are  higher  in  both  analyses.  Hence  the 
ft  versus  crack  Length  curves  for  various  applied  loads  will  have  the  same 
shape  as  shown  in  Figure  9 for  the  intact  stringer  case,  however*  the  /?  values 
will  be  higher  at  the  same  crack  length  and  applied  stress*  Thus  in  the  case 
of  the  broken  stringer  the  first  point  of  instability  will  be  at  some  stress 
lower  t nan  p/F  t = 0*545  and  crack  arrest  will  occur*  The  tested  crack 
lengths  for  broken  and  intact  stringer  are  similar  and  hence  values  at 

failure  for  these  two  panels  would  be  expected  to  be  similar*  Using  the 
predicted  value  p/F ty  - 0*645  and  the  elastic-plastic  analysis  for  a half 
crack  length  of  2*75" Inches  (Figure  6,  intact  stringer)  the  /I  value  at 


Using  this  /J  value  and  the  elastic-plastic 


failure  is  31,7 


analysis  for  the  broken  stringer  panel  with  a half  crack  length  of  2*75 
inches  (Figure  6)  the  predicted  failure  stress  is  p/Fty  - 0.553.  The 

actual  failure  load  of  the  -23  panel  was  107*2  kips  whereas  the  predicted 
failure  Load  corresponding  to  p/F  ty  = 0.558  is  106,74  kips.  Thus  the 

predicted  load  is  within  0.44  percent  of  the  actual  failure  load, 

2*3  TITANIUM  PANEL 

To  determine  the  influence  of  material  and  toughness  (through  variation 
of  resistance  curve)  on  the  residual  strength  of  fatigue  cracked  panels*  an 
all  titanium  panel  was  tested  and  residual  strength  determined*  The  skin 
materiaL  consisted  of  a beta  mill  annealed  Ti”6Al-4V  and  formed  angles  oi 
annealed  Ti - 6AL - 6V-2Sn.  The  construction  once  again  as  with  the  previous 
panels  represents  a typical  lower  surface  wing  construction  sized  for  a 
fighter  type  aircraft. 

2,3*1  Description  of  Test  Panel 

Figure  L indicates  the  drawing  of  the  titanium  test  panel  and  is 
identified  as  a -21  assembly,  however,  the  term  titanium  panel  is  used  in 
the  report*  Overall  skin  width  was  38*5  inches  with  a panel  Length  of 
approximately  90  inches.  Spacing  between  rivet  and  angle  centerlines  was 
5*5  inches*  This  spacing  was  identical  to  that  of  the  previous  wing  panels 
and  results  in  a six  bay  panel  symmetric  about  the  panel  centerline. 
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The  titanium  angles  were  fabricated  from  0.064  inch  thick,  Ti - 6A1 -6V -2Sn 
(mill  annealed)  by  passing  through  a heated  die.  The  material  is  the  same  as 
that  previously  characterized  in  Reference  2 and  no  change  in  material 
properties  of  the  attached  leg  were  noted  due  to  the  hot  forming  operation. 
Dimensions  of  the  angle  are  shown  in  Figure  1,  (Section  B-B)  . As  with  the 
previous  aluminum  wing  panels  the  angles  were  attached  through  the  one  inch 
leg  by  3/ 16-inch  diameter,  flush  head  steel  HI-LOK's  (19PB-3)  at  1.3  inch 
pitch  in  the  test  area.  Listed  in  Table  V are  the  panel,  skin,  and  angle 
cross  sectional  thicknesses,  areas,  and  material  identification.  It  will  be 
noted  that  the  stiffening  ratio  of  this  panel  is  within  0.2  percent  of  the 
-15  aluminum  wing  panel  (see,  e,g,s  Table  1), 


TABLE  V TITANIUM  SKIN -ANGLE  STIFFENED  WING  PANEL  AREAS  AND  MATERIAL 
IDENTIFICATION 


MEASURED  THICKNESS 

( Inches) 

CROSS  SECTIONAL  AREAS 

2 

(Inches  ) 

SKIN 

MATERIAL/ 

AND  L.D. 

(See  Ref.  2) 

PERCENTAGE 

OF  TOTAL 

ANGLE  AREA 

TO  SKIN 

AREA 

Angle 

Skin 

Total  Angies 

Skin 

.064 

.058 

.6524 

2.233 

Beta  Mill 
T1-6AL-4V/ 

TO  #3 

29.27. 

Total  Panel  Cr 
Sectional  Area 
2.8854  in2 

OSS 

Indicated  in  Figure  l is  a starter  slot  5-00  inches  in  overall  length. 

As  with  the  previous  panels  20  strain  gages  were  positioned  as  shown  in 
Figure  27.  In  addition,  the  beam  clip  gage  was  mounted  on  the  skin  side 
to  record  crack  opening  displacement  at  the  crack  and  panel  centerline. 

2.3. 2 Finite  Element  Modeling  and  Analysis 

The  analyses  of  the  all  titanium  panel  was  performed  in  a manner  similar 
to  the  intact  stringer,  all  aluminum  panel  (-15).  Both  elastic  and  elastic- 
plastic  analyses  were  performed  for  various  crack  lengths.  The  finite  element 
modeling,  elastic,  and  elastic-plastic  analyses  of  the  panels  are  discussed  in 
the  following  subsections. 


41 


w 


e 

«d 

EXi 

E 

D 

a 

rs 

j-i 

■*rJ 

H 


(13 

C 

O 


d 

u 

o 

-I 

OJ 

CkD 

<d 

a 


cd 

lul 


cn 


r-' 

cm 

to 

u 

3 

&d 

u. 


^ 1 L 1 1 ll  1 


42 


2 b 3 , 2 . 1 Finite  Element  Modeling  of  the  All  Titanium  Panel 


The  finite  element  model  used  to  analyze  the  aluminum  panel  C-15)  was 
adopted  to  treat  the  titanium  panel.  The  size  of  the  angles  and  skin  thickness 
for  this  panel  were  different  from  those  of  the  aluminum  panel.  The  lx,  y) 
co-ordinates  of  all  grid  points  in  the  finite  element  model  were  kept  the  same 
as  for  the  aluminum  panel  and  only  the  z (dimension  perpendicular  to  the  plane 
of  the  panel,  see  Figure'  2)  coordinates  required  changing.  Changes  were 
made  in  the  thicknesses  and  areas  of  the  triangular  membrane  and  rod  elements 
of  the  finite  element  model  to  comply  with  the  dimensions  used  for  the  titanium 

pane  I . 

2.3,2. 2 Elastic  and  Elastic-Plastic  Analyses  of  All  Titanium  Panel 

An  elastic  analysis  of  the  titanium  panel  with  intact  central  stringer 
was  performed  for  a half  crack  length  of  2.875  inches  The  contours  shown  n 
Figure  3 were  used  to  evaluate  J.  The  variation  of  /J  with  p/Fty  is  shown  in 

Figure  28* 


The  elastic-plastic  analysis  was  conducted  assuming  a Dugdale  type  strip 
plastic  zone  ahead  of  the  crack  tip.  The  variation  of  /J  with  p/Fty*  using 
Dugdale  type  elastic -plastic  analysis,  is  shown  in  Figure  28  for  a half  crack 
length  of  2.875  inches.  At  an  applied  stress  of  80  percent  of  yield  the 
values  of  /J  given  by  an  elastic-plastic  analysis  are  18  percent  higher  than 
those  given  by  an  elastic  analysis.  To  obtain  ,/J  versus  crack  length  curves 
for  various  applied  loads  required  for  the  residual  strength  prediction,  five 
half  crack  lengths  were  selected:  1.75,  2.875,  4.75,  5.5,  and  6.0  inches. 

The  plot  of  /J  versus  applied  stress  is  shown  in  Figure  29  for  these  physical 
crack  lengths.  These  /J  values  have  been  cross  plotted  in  Figure  30  to  show 
the  variation  of  /J  with  crack  Length  for  various  applied  loads  (stress). 

From  Figure  30  it  can  be  seen  that  the  /J  values  increase  with  crack  length 
to  a half  crack  length  of  approximately  4.73  inches  and  then  decrease  wiLh 
additional  crack  length.  The  /J  values  reach  a minimum  value  at  the  stringer 
(angle)  Located  at  5.5  inches  from  the  crack  and  panel  centerline.  When  the 
crack  is  beyond  the  centerline  of  the  stringer,  the  values  increase  again 
with  crack  length. 


2.3.3  Experimental  Results  and  Discussion 

In  References  1 and  2 the  data  necessary  to  evaluate  crack  growth 


resistance  based  on  Kr  versus  da  were  presen 


ted  for  the  T1-6AL-4V  Beta  Mill 


Annealed  material.  However,  those  data  were  not  presented  based  on  the 
proposed  failure  criterion,  i.e.,  /Jg  versus  ua.  As  a prelude 
tical  and  experimental  evaluation  of  the  residual  strength  of  this  titanium 
panel  these  data  will  be  presented  in  the  next  section.  References  and 
should  be  consulted  for  details  of  how  these  data  were  obtained  and  Jhe 
associated  heat  treatment,  microstructure  and  mechanical  property  data  o 
this  material- 


Square  Root  of  Jf  /J  /in,  - lbs 


Figure  28.  Elastic  and  Elastic  Plastic  Analysis  of  All  Titanium  Panel 
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Figure  29,  Square  Root  of  J as  a Function  of  Applied  Stress  for  Varying  Crack  Lengths  - 
All  Titanium  PaneL 
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Square  Root  of  J,  /J  ( ilb — I — 


F i gu  re  30 , 


Square  Root  of  J as  a Function  of  Crack  Length  - Titanium  Panel 
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2, 3. 3*1  /Jr  Crack  Growth  Resistance  Data  for  Beta  Mill  Annealed  Ti-6AL-4V 


Two  thicknesses  (*053  and  *187  inches)  and  two  fracture  plane  orienta- 
tions ( LT  and  TL)  were  examined,  using  the  crack  Line  wedge  loaded  (CLWL) 
specimen  geometry*  Figure  31  indicates  the  */Jr  versus  t± ap^y  _ values  for 
the  LT  direction  and  Figure  32  for  the  TL  direction  in  the  thinner  gage  and 
Figures  33  and  34  for  the  thicker  gage  material  in  the  LT  and  TL  direction 
respectively*  Consistency  of  data  is  quite  evident  for  the  thinner  gages* 

The  thicker  gage  material  in  the  LT  direction  shows  evidence  of  banding  which 
is  reflected  from  crack  deviation  from  a plane  normal  to  the  loading  axis 
for  this  material  (see  Reference  2 for  specific  discussion  of  this  behavior), 

2 , 3 * 3 , 2 Fatigue  F recracking,  Strain  Data  and  Fracture  of  Titanium  Panel 

Ihe  titanium  panel  was  precracked  In  the  same  manner  as  the  previous 
aluminum  panels.  Sinusoidal,  constant  amplitude,  tension- tension  loading  was 
employed  and  an  ini tial  maximum  stress  intensity  of  approximately  20  ksi 
/inch  and  then  25  ksi  /inch  selected  for  fatigue  loading*  Subsequent  loading 
was  at  lower  levels  to  sharpen  the  crack*  Table  VI  summarizes  the 

fatigue  loading  and  total  crack  length  for  this  titanium  panel.  The  final 
fatigue  crack  lengths  were  measured  after  fracture. 


TABLE  VI  FATIGUE  PRECRACKING  AND  CRACK  LENGTH 
DATA  - TITANIUM  PANEL 


FATIGUE 

LOADS 

NUMBER  OF 

CYCLES 

FATIGUE 

. CRACK  LENGTH(S) 

MAXIMUM 

(kips) 

MINIMUM 

(kips) 

LEFT 
( Inch) 

RIGHT 

(Inch) 

TOTAL 

(Inch) 

20.0 

3.0 

3,000 

0 

0 

0 

33.5 

u.o 

1 8 , 200 

* 11 

.08 

5. 19 

27.0 

3.0 

12,050 

*35 

.31 

5. 66 

Figures  35  and  36  show  the  front  (skin)  side  view  and  angle  side  view 
of  the  all  titanium  panel* 

After  fatigue  cycling  and  zero  balancing  of  all  strain  gages  and  displace- 
ment gage  a strain/ Load  survey  was  taken  by  loading  to  specific  load  intervals* 
These  data  are  noted  in  Table  VII*  Figure  27  has  been  repeated  as  Figure  37 
following  the  strain  table  for  convenience  in  locating  gage  positions* 
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Figure  31.  /J^  Crack  Growth  Resistance  Data  - T1-6A1-4V 
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Figure  35,  Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened,  Lower 
(All  Titanium  Panel) 
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Figure  36.  Overall  View  of  Angle  Stiffened  Side  of  Wing  Panel  Angle  Stiffened, 
Lower  (All  Titanium  Panel) 
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TABLE  VII  STRAIN  GAGE  DATA,  TITANIUM  WING  PANEL  STIFFENED , LOWER 
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GAGES  ON  UPSTANDLNG  LEG  OF  ANGLE 


Load  versus  total  displacement  gage  output  is  shown  in  Figure  38*  Four 
distinct  loadings  are  noted*  Strain  data  were  recorded  for  the  last  three 
load  excursions,  the  first  of  these  involved  a strain  survey  from  zero  load  to 
80  kips . After  unloading, loading  was  re-initiated  using  ramp  loading  of 
100,000  pounds  per  minute  until  some  visible  crack  extension  took  place. 

This  extension  0.10  inch)  is  noted  on  both  the  strain  (Table  VI t)  and 
displacement  data  (Figure  38.)  at  a load  of  142,4  kips.  At  this  point  the 
load  was  again  reduced  to  zero  and  ramp  loading  once  again  initiated  after 
measurement  of  the  crack  extension. 

During  the  loading  to  failure  the  crack  extended  slowly  to  the  adjacent 
stringers  at  which  time  an  attempt  was  made  to  hold  the  load.  However,  the 
arrest  was  of  such  short  duration  (<  3 seconds)  compared  to  the  loading  rate 
(100,000  pounds/minute)  that  the  instruction  to  the  operator  to  hold  the 
load  was  futile.  Final  fracture  for  this  titanium  panel  was  at  220,6  kips. 

Figures  39  and  40  show  the  fracture  surface  of  the  titanium  panel  from 
the  skin  and  angle  sides.  Due  to  the  large  grained  structure  for  the  skin 
material  the  zig-zag  appearance  was  not  unexpected.  This  behavior  for  the 
Beta  Mill  annealed  condition  was  quite  common  to  this  material  as  noted  in 
Reference  2.  Fracture  of  the  angles  essentially  occurred  along  the  line  of 
center  fasteners  with  one  exception  as  shown  in  Figure  40. 

2 * 3 , i , 3 Comparison  of  Experimental  and  Analytical  Strains 

Figure  37  shows  the  location  of  strain  gages  on  the  titanium  panel. 
Tabulated  in  Table  VII  are  the  strain  gage  readings  for  various  applied  loads 
at  a crack  Length  of  5,66  inches.  Figure  41  shows  the  experimental  and 
analytical  variation  of  strain  in  the  central  stringer  with  increasing 
distance  from  the  crack.  Good  correlation  is  noted  between  experimental 
strains  and  those  obtained  from  the  elastic  analysis.  The  experimental  data 
of  Figure  41  is  cross  plotted  in  Figure  42  to  show  the  variation  of  strain 
with  applied  stress  for  this  panel.  Good  correlation  is  obtained  between 
elastic  analysis  and  experimental  data  to  an  applied  stress  of  28  ksi*  At 
higher  applied  stresses  the  experimental  strain  in  the  stringer  close  to  the 
crack  plane  is  lower  than  that  given  by  an  elastic  analysis.  Away  from  the 
crack  plane,  the  correlation  between  experimental  strains  and  those  obtained 
by  elastic  analysis  is  good  to  an  applied  stress  of  60  ksi, 

2,3,4  Residual  Strength  Prediction 

The  residual  strength  prediction  procedure  discussed  in  Section  2.2,4 
(see  also  Reference  1)  was  used  to  predict  the  residual  strength  of  the  aLi 
titanium  panel.  The  analytical  stresses  in  the  central  stringer  (angle) 
with  a half  crack  Length  of  2,875  inches  are  shown  in  Figure  43,  The 
ultimate  strength  of  the  stringer  material  (Ti - 6A 1 - 6V - 2Sn ) was  determined 
through  test  to  be  167.3  ksi.  Using  this  ultimate  strength  value  and  the 
Dugdale  type  el ast ic - plastic  analysis  stringer  stresses  were  determined  as 
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Total  Displacement f 6 (inches) 


Figure  39,  View  of  Fracture  Surface  From  Skin  Side,  Wing  Panel  Angle 
Stiffened,  Lower  t A 1 1 Titanium  Panel) 
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Figure  40.  View  of  Fracture  Surface  From  Stringer  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (All  Titanium  Panel) 
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Figure  41.  Strain  in  Central  Stringer  for  All  Titanium  Panel,  Angle  Stiffened,  Lower 
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Strain  In  Central  Stringer,  (inches/inch 


Figure  42.  Strain  in  Central  Attached  Angle  Away  From  the  Crack  Plane 
as  a Function  of  Applied  Stress  for  All  Titanium  Panel 
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43.  Elastic  and  Elastic-Plastic  Stress  in  Central  Stringer 
as  a Function  of  Applied  Stress  - All  Titanium  Panel 
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shown  in  Figure  43.  The  failure  of  the  central  stringer  tangle)  is  predicted 
at  p/Fty  - 0.735.  Now  consider  the  /j  versus  crack  length  plot  for  various 

appLied  stresses  as  shown  in  Figure  44  for  the  all  titanium  panel.  Super- 
imposed on  this  plot  is  the  /J  g versus  A ap^y  resistance  curve  of  the  Beta 

Mill  annealed  Ti tanium-6Al-4V  skin  material  plotted  at  a physical  half  crack 
length  of  2.83  inches.  It  can  be  noted  from  this  Figure  at  the  stringer 
critical  stress  of  p/F^y  = 0.735  (determined  from  Figure  43)  there  would  be 

considerable  slow  tear  in  the  skin.  Therefore  the  stringer  stresses  shown  in 
Figure  43  will  not  be  valid  at  failure  as  the  crack  length  would  have  increased 
considerably  prior  to  stringer  failure.  From  Figure  44*  the  resistance  curve 
of  the  material  plotted  at  the  initial  fatigue  half  crack  length  of  2.83 
inches  (initial  half  crack  length  in  the  tested  panel)*  is  tangential  to  the 
/j  curve  of  the  panel  at  an  applied  stress  of  p/Fty  - 0,66.  Thus  the  first 

point  of  crack  instability  (rapid  crack  extension)  occurs  after  an  initial 
period  of  slow  crack  growth.  However,  from  Figure  44  it  can  be  noted  that 
immediately  after  the  point  of  instability  the  curve  for  the  panel  drops 
lower  than  the  resistance  curve  of  the  material  (L.e.,  beyond  this  point  the 
resistance  of  the  material  is  higher  than  /j  developed  in  the  stressed  panel) , 
Hence  the  running  crack  will  be  arrested.  The  crack  will  become  arrested 
where  the  next  stringer  (angle)  is  connected  to  the  skin.  The  resistance 
curve  of  the  material  is  now  replotted  at  a half  crack  length  of  5.5  inches 
(distance  from  centerline  of  the  stringer  to  panel  centerline)  where  the  crack 
became  arrested.  For  this  crack  length  the  AJ  curve  of  the  paneL  becomes 
tangent  to  the  resistance  curve  at  an  applied  stress  of  p/F^y  - 0,69.  Hence, 

the  second  point  of  crack  instability  occurs  at  this  stress  and  the  crack  runs 
catastrophically.  Beyond  this  point  of  instability  upon  further  Loading  the 
/J  values  in  the  panel  are  much  higher  than  the  resistance  curve  of  the  skin 
material  and  therefore  no  possibility  exists  for  crack  arrest.  It  may  be 
noted  that  the  difference  in  the  stress  at  the  first  point  of  instability 
(i.e,,  P^ty  = 0.66)  and  at  the  second  point  of  instability  (i.e.,  p/FL^,  - 

0,69)  is  small.  Hence,  under  increasing  applied  load,  the  crack  will  only 
momentarily  arrest. 


The  titanium  panel  analyzed  above  was  tested  to  failure.  The  slow  tear 


in  this  panel  started  at  a load  of  142.4  kips 


(*1r0-43)' 


At  this  load 


the  slow  tear  at  each  end  was  0.1  inch  (see  Table  VII),  The  first  point  of 
instability  (rapid  crack  extension)  was  at  a load  Lower  than  220.6  kips 

The  crack  was  arrested  at  the  rivet  hole  in  the  adjoining 


stringer(s)  for  only  a short  duration  and  before  the  increasing  applied  load 


could  be  held  constant  panel  failure  occurred.  The  time  lag  between  the 
first  and  second  points  of  instability  was  less  than  five  seconds.  For  this 
panel  the  predicted  first  point  of  instability  is  at  a load  of  217.6?  kips 

= 0.66]  and  final  failure  of  the  panel  occurred  at  a load  of  227,56  kips 

lFty  ) 

= 0*69^J  . This  load  is  3.15  percent  higher  than  the  measured  failure 
load  of  220.6  kips* 
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Figure  44.  Residual  Strength  Prediction  Plot  for  All  Titanium  Panel 
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2.4  INTACT  STRINGER  - THICK  SKIN 


The  Intent  of  testing  this  "thick  skin"  panel  was  to  determine  the  effect 
of  mixed  mode  fracture  behavior  on  residual  strength  prediction.  The  material 
and  thickness  selected  for  this  panel  was  chosen  so  that  a mixed  mode  slow 
tear  behavior  would  result.  Due  to  the  increased  skin  gage  a rather  low 
percentage  stiffening  occurred  but  the  overall  objective  of  a mixed  mode 
slow  tear  was  accomplished. 

2 , 4 , L Description  of  Test  Panel 

Figure  1 gives  the  detailed  geometry  of  this  "thick  skinned"  test  panel. 
It  is  designated  as  a -1  Assembly  with  a test  section  width  of  38,5  inches. 

The  7075-T73  skin  material  had  symmetrical  angle  sections  (AND1G133 - 1002 , 
7075-T6  extrusions)  attached  through  the  leg  by  3/16-inch  diameter,  flush 
head  steel  HI-LOKS  at  1.3 -inch  pitch.  A bay  width  of  5.5  inches  resulted 
in  a six  bay  panel  symmetric  about  the  panel  centerline, 

TabLe  VX11  lists  the  panel,  skin,  and  angle  c ross - sec t iona L thickness, 
areas  and  material  identification.  At  a 17  percent  stiffening  ratio  this 
panel  would  be  considered  to  be  a Lightly  stiffened  panel. 


TABLE  VIII  THICK  SKIN  (-1)  - ANGLE  STIFFENED  WING  PANEL 
AREAS  AND  MATERIAL  IDENTIFICATION 


MEASURED  THICKNESS 
(Inches ) 

CROSS  SECTIONAL  AREAS 

( Inches^ ) 

SKIN  MATERIAL/ 
AND  I.D. 

(See  Ref.  2) 

PERCENTAGE 

OF  TOTAL 

ANGLE  AREA 

TO  SKIN  AREA 

Angle 

Skin 

Total  Angles 

Skin 

. 094 

. 193 

I .281 

7.4305 

7075-T73/ 

5-2 

17.2 

Total  Panel  Cross 
Sectional 

Area  - 8.7115  in , ^ 

The  initial  starter  slot  half  length 
of  8,50  inches.  Twenty  strain  gages  were 
As  with  the  other  wing  panels  a beam  type 
side  to  record  crack  opening  displacement 


was  4,25  inches  or  an  overall  length 
positioned  as  shown  in  Figure  45. 
clip  gage  was  mounted  on  the  skin 
at  the  crack  and  panel  center  Line, 
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GAGES  ON  ANGLE  SIDE  ONLY 
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Figure  45*  Strain  Gage  Locations  - Thick  Skin 


2.4.2  Finite  Element  Modeling  and  Analysis 


The  finite  element  model  employed  to  analyze  the  -15  panel  was  adopted 
for  use  in  this  thick  skin  panel.  Changes  were  required  in  both  skin  thickness 
and  stiffener  areas  to  correspond  to  the  actual  skin  thickness  and  stiffener 
areas  of  this  (-1)  panel.  The  lx,  y)  co-ordinates  of  all  grid  points  in  the 
finite  element  model  of  the  -l  panel  were  the  same  as  previously  employed 
for  the  -15  panel.  The  z co-ordinates  of  grid  points  connecting  rod  and 
membrane  elements  l representing  stiffeners)  were  changed  to  correspond  to  the 
skin  thickness  and  stiffener  dimensions  of  the  -1  panel. 

An  elastic  analysis  was  conducted  at  a half  crack  length  of  4.5  inches. 

The  contours  shown  in  Figure  3 were  used  to  evaluate  J.  The  variation  of 
with  p/Fty  is  shown  in  Figure  46. 

The  elastic- plastic  analysis  of  this  panel  was  performed  assuming  a 
Dugdale  type  strip  plastic  zone  ahead  of  the  crack  tip.  The  variation  of  ✓ J 
with  p/Fty  is  shown  in  Figure  46  for  a half  crack  length  of  4.5  inches.  At  an 
applied  stress  of  80  percent  of  skin  material  yield,  the  fJ  value  given  by  an 
elastic  analysis  is  7 percent  lower  than  that  given  by  assuming  a Dugdale 
model  analysis.  An  elastic-plastic  analysis  for  various  assumed  crack  lengths 
was  npt  performed  on  this  panel.  The  stresses  in  the  central  stringer  of  the 
panel  for  elastic  and  elastic -plastic  analyses  are  shown  in  Figure  4?. 

2.4.3  Experimental  Results  and  Discussion 

This  -l  Assembly  was  the  first  panel  fabricated  and  tested  in  the  wing 
panel  configuration.  During  one  of  the  unloading  sequences  the  specimen  was 
inadvertently  placed  into  compression  loading  which  was  sufficient  to  produce 
a permanent  set  in  the  angle  sections  in  the  test  area.  This  resulted  in  a 
situation  which  adversely  affected  the  resulting  residual  strength  prediction. 
However,  the  strain  and  fracture  data  will  be  presented  here  with  additional 
discussion  to  follow. 

2.4.3. 1 Fatigue  Precracking,  Strain  Data  and  Fracture  of  Thick  Skin  Panel 

Fatigue  pre cracking  of  the  -1  Assembly  was  accomplished  in  a manner 
similar  to  the  previous  wing  panels.  Sinusoidal,  constant  amplitude,  tension- 
tension  loading  was  employed  at  a maximum  stress  intensity  level  of  15  ksi 
/inch.  Table  IX  summarizes  the  fatigue  loading  and  final  crack  lengths  for 

this  thick,  7075-T73  panel. 

Figures  48  and  49  indicate  skin  and  angle  side  views  of  the  -l  Assembly 
in  the  500  kip  test  machine  prior  to  testing. 

After  fatigue  cycling  and  zero  balancing  of  all  strain  gages  and  dis- 
placement gage  a strain/ load  survey  was  taken  by  loading  to  specific  load 
intervals.  Gages  1,  3,  and  6 were  monitored  on  tape  as  well  as  by  direct 
strain  readings  to  compare  strain  data  for  the  two  systems.  These  data  are 
given  in  Table  X.  It  can  be  seen  from  these  data  that  as  with  the  comparison  of 
strain  data  for  the  two  strain  recording  systems  in  Reference  1 the  correlation  be- 
tween the  two  is  quite  good.  After  the  tear  to  the  adjacent  angles  the  panel  was 
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Figure  46.  Elastic  and  Elastic-Plastic  fZ  Values  for  Thick  Skin, 
-1  Panel  as  a Function  of  Applied  Stress 
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Figure  U7.  Elastic  and  Elastic-Plastic  Stress  in  Central  Stringer  as 
Function  of  Applied  Stress  -1  Panel 
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inadvertently  placed  in  compression.  The  magnitude  of  the  permanent  set  placed 
in  the  panel  can  be  seen  from  the  strain  data  at  zero  load  (after  tear  to 
angles  - Table  X)  where  the  strain  at  gage  number  10  l s quite  large,  Gage  10 
is  located  on  the  upstanding  leg  of  the  angle  near  the  crack,  see  e,g,  , Figure 
45,  Starting  from  zero  load  the  initial  strain  survey  went  from  zero  to  140 
kips* 


TABLE  IX  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA  * 
THICK  SKIN  PANEL 


FATIGUE  LOADS 

NUMBER  OF  CYCLES 

FATIGUE  CRACK  LENGTH(S) 

MAXIMUM 

(kips) 

MINIMUM 

(kips) 

LEFT 

(Inch) 

RIGHT 

(Inch) 

TOTAL 

(Inch) 

44.0 

4.4 

11,350 

4.34 

4.50 

8.84 

Load  versus  total  displacement  at  the  crack  centerline  are  shown  in 
Figure  50  for  the  load  survey  and  slow  tear  and  the  final  fracture  in 
Figure  51.  It  can  be  seen  that  for  the  loading  survey  to  140  kips  the  zero 
on  the  deflection  scale  was  stepped  to  aid  in  data  reduction.  On  the  loading 
segment  which  resulted  in  tear  to  the  adjacent  angles  (see  third  curve  of 
Figure  50,  two  significant  steps  can  be  noted.  They  resemble  the  pop- in 
behavior  of  unreinforced  panels  of  7075-T73  aluminum  of  this  thickness. 

After  crack  arrest  the  load  was  held  at  251,1  kips,  strain  readings  recorded 
and  the  panel  unloaded  to  zero.  At  this  point  the  panel  was  inadvertently 
placed  into  compression  at  which  point  noticeable  "straining"  of  the  angles 
in  the  panel  centerline  occurred.  The  permanent  set  which  occurred  during 
this  sequence  is  quite  apparent  from  the  initial  part  of  the  load -d i splace - 
ment  trace  to  failure  in  Figure  51,  Tne  reverse  curvature  is  common  to  center 
cracked  panels  which  have  undergone  prior  inplane  buckling  neat  the  crack. 

Upon  further  loading  fracture  occurred  from  this  -l  Assembly  at  337.8  kips. 

Figure  52  shows  the  angle  side  of  the  panel  after  fracture.  In  Figure  53 
the  close  up  of  the  fracture  surface  indicates  that  with  the  exception  of  two 
of  the  seven  angles  fracture  occurred  at  the  first  line  of  fastener  away  from 
the  crack  line.  In  addition  the  permanent  set  due  to  the  compression  load 
is  quite  evident  in  the  central  and  adjacent  angles.  The  broken  piece  of 
angle  indicates  the  severity  of  the  permanent  set  which  occurred  in  this 
panel.  It  can  also  be  seen  that  during  the  tearing  to  the  angles  both  mixed 
mode  and  shear  crack  progression  occurred  until  the  crack  arrested  at  the 
fastener  hole(s).  Upon  additional  loading  to  fracture  shear  crack  progression 
predominates . 
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Figure  48.  Overall  View  of  Skin  Side  of  Wing  Panel  Angle  Stiffened,  Lower 
(Thick  Skin) 
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Figure  49.  Overall  View  of  Angle  Stiffened  Side  of  Wing  Panel  Angle  Stiffened, 
Lower  (Thick  Skin) 
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TABLE  X STRAIN  GAGE  DATA,  -1  THICK  SKIN  WING  PANEL  ANGLE  STIFFENED,  LOWER 
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250  “ •*  1 f ~~ j— j — S_L  Holding  Load  After  Tear  to  Angles-: 
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Total  Displacement  5 (inches) 


-i  r in.  i 

Figure  52.  View  of  Panel  Fracture  From  Stringer  Side,  Wing  Panel  Angle 
Stiffened,  Lower  (Thick  Skin) 
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Figure  53.  Close-up  View  of  Fractured  Angles,  Wing  Panel  Angle  Stiffened, 
Lower  (.Thick  Skin) 
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2 . 4 . 3 , 2 Comparison  of  Experimental  and  Analytical  Strains 


Figure  54  shows  the  location  of  strain  gages  on  the  -l  panel.  Tabulated 
in  Table  X are  the  strain  gage  readings  for  various  applied  loads  at  a fatigue 
crack  length  of  8.84  inches.  Figure  55  shows  the  experimental  and  analytical 
variation  of  strain  in  the  central  stringer  with  increasing  distance  up  the 
central  angle  from  the  crack  at  a half  crack  length  of  4.42  inches.  Good 
correlation  is  noted  between  experimental  data  and  the  elastic,  finite  element 
analysis.  The  experimental  data  of  Figure  55  has  been  cross  plotted  in 
Figure  56  to  show  the  variation  of  strain  with  applied  stress.  Good  correla- 
tion is  obtained  between  the  experimental  strain  data  and  the  elastic  analysis 
to  an  applied  stress  of  approximately  18  ksi.  For  higher  stresses  the  strains 
in  the  central  stringer  near  the  plane  of  the  crack  are  lower  than  the  elastic 
strains.  Similar  behavior  was  observed  for  the  titanium  panel, 

2.4.4  Residual  Strength  Prediction 

The  elastic -plastic  analysis  of  this  panel  was  not  performed  for  various 
crack  lengths.  Hence,  curves  between  /J  and  several  physical  crack  lengths 
for  various  applied  loads  could  not  be  pLotted.  It  may  be  noted  that  the 
structural  arrangement  of  the  thick  skin,  -1  panel  and  the  thin  skin  -15 
panel  are  the  same.  Thus  a plot  of  the  trend  in  /j  versus  a for  the  -1 
panel  will  be  similar  to  that  for  the  -15  panel  (Figure  9),  however,  in 
this  case  the  /J  values  will  be  different  for  the  same  applied  stress  and 
crack  length.  For  a half  crack  length  of  4,45  inches  the  /j  values  for  the 
-1  panel  were  determined  to  be  between  17  to  19  percent  higher  than  those  for 
the  -15  panel.  Thus  for  the  residual  strength  prediction  of  the  - 1 panel  it 
was  assumed  that  the  /J  values,  for  the  various  crack  lengths  and  applied 
stress,  were  18  percent  higher  than  those  for  the  -15  panel.  Based  on  this 
assumption,  /T  versus  a curves  for  the  -1  panel  have  been  plotted  in  Figure 
57  using  an  18  percent  increase  to  the  curves  of  Figure  9.  These  curves 
were  then  used  to  determine  the  skin  critical  behavior  in  the  panel* 

The  stresses  in  the  central  stringer  for  the  -l  panel  are  shown  in 
Figure  47,  The  ultimate  strength  of  the  stringer  material  (7073-T6 
extrusions)  was  determined  by  test  to  be  87.2  ksi.  Using  this  ultimate 
stress  value  and  the  Dugdale  type  elastic  - plastic  analysis  shown  in  Figure 
47,  the  failure  of  the  central  stringer  angle  is  predicted  at  an  applied 
stress  ratio  of  p/Fty  = 0.455.  The  /J  versus  crack  length  plots  for  various 
applied  stresses  are  shown  in  Figure  57.  Superimposed  on  this  plot  is  the 
/Jr  versus  aap^.  resistance  curve  of  the  0,195  inch  thick  7075-773  (LX)  skin 

material  plotted  at  a physical  half  crack  length  of  4,42  inches.  it  can  be 
seen  that  at  the  stringer  critical  stress  p/Ft  - 0.455,  there  will  be  slow 

tear  in  the  panel  equal  to  approximately  0,8  inch  at  each  end  of  the  crack. 
Hence  the  stringer  stresses  obtained  by  a Dugdale  analysis  for  a half  crack 
length  of  4.42  inches  will  not  be  valid  for  this  panel.  From  Figure  57 
it  is  seen  that  at  an  applied  stress  ratio  p/Fty  - 0,5  the  crack  will  slow 

tear  to  the  rivet  hole  connecting  the  adjacent  stiffener  to  skin,  or  a point 
very  ciose  to  the  angle  centerline.  It  may  be  noted  that  no  instability  of 
the  crack  occurs  before  the  crack  slow  tears  to  the  hole  due  to  the  /J  curves 
falling  below  the  /J  ^ curve.  The  /Jr  curve  was  rep  lot  ted  at  the  half  crack 
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Strain  in  Central  Stringer,  Unches/inch  x 10^) 


Figure  55-  Strain  in  Central  Stringer  for  Thick  Skin,  -i  Panel, 
Angle  Stiffened,  Lower 
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Figure  56.  Strain  in  Central  Stringer  as  a Function  of  Applied  Stress  -l  Panel 
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length  of  5_5  inches  or  at  the  centerline  of  the  stringer.  For  this  crack 
Length  y/J  curve  of  the  panel  becomes  tangential  to  the  resistance  curve  at 
an  applied  stress  ratio  of  p/Fty  = 0.55,  Hence  the  instability  should  occur 
at  this  stress  and  the  crack  will  start  running  catast rophical Ly- 

The  panel  analyzed  above  was  tested  to  failure.  Upon  loading  the  crack 
tore  to  the  rivet  hotels)  at  a Load  of  241  kips  or  p/F^y  — 0,411  (predicted 

p/F  r - 0,50) - After  the  slow  tear  to  the  rivet  hole,  the  panel  was  inad- 

verten L 1 y placed  into  compression  loading  which  was  sufficient  to  produce  a 
permanent  set  in  the  angle  sections  in  the  test  area.  Due  to  this  action  the 

central  stringer  was  able  to  Lake  a much  higher  load.  Ihe  final  failure  of 

the  panel  occurred  at  an  applied  load  of  337.8  kips.  Ihis  corresponds  to  a 
predicted  failure  load  of  (p/Fty  = 0.55)  322. A56  kips.  The  predicted  failure- 
load is  within  A. 7 percent  of  the  measured  fai lure- load . However,  it  must  be 
noted  that  the  /5  versus  a curves  for  this  panel  shown  in  Figure  57  were  not 
obtained  by  rigorous  analysis  and  therefore  represent  approximate  values. 

Thus,  this  predicted  failure-load  can  be  only  approximate.  It  should  be 
noted  that  the  failure  of  this  paneL  cannot  be  analyzed  by  the  proposed 
technique  due  to  the  compression  buckling  which  occurred  in  the  test  machine. 
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Ill  STIFFENED  BIAXIALLY  LOADED  PANELS 


In  Section  II  the  case  of  uniaxial Ly  Loaded  panels  were  examined  which 
would  typically  represent  sections  of  lower  aircraft  wings.  It  is  known  that 
in  most  fuselage  structure  more  complex  loading  conditions  predominate  and 
muL t i -axial  loading  prevails.  In  light  of  these  anticipated  Loading  conditions 
it  was  believed  necessary  to  simulate  and  demonstrate  the  prediction  technique 
for  some  form  of  multi-axial  loading.  The  first  decision  involved  choice  of 
loading  and  the  second  the  representative  structure. 

Choice  of  loading  was  biaxial  and  the  structure  was  representative  of 
the  upper  fuselage  of  a large  transport/bomber  type  aircraft.  The  number  of 
panels  tested  (two)  permitted  only  a cursory  examination  of  the  problem  of 
biaxial  loading  on  residual  strength.  However,  trends  were  observed  and  some 
tentative  evaluation  of  the  proposed  anaLysis  was  possible.  This  evaluation 
will  be  discussed  in  subsequent  subsections. 

3.1  GENERAL  INTRODUCTION 

Ihe  influence  of  biaxial  loading  on  the  crack  propagation  and  fracture 
strength  of  a cracked  structure  has  attracted  considerable  attention  in  recent 
years.  The  biaxial  Loading  represents  a realistic  state  of  stress  in  a 
structure.  The  majority  of  stress  analysis  of  cracked  structures  available 
in  the  literature  deal  with  uniaxial ly  loaded  cases  with  a loading  direction 
perpendicular  to  the  plane  of  the  crack.  Research  has  been  limited  in  the 
area  of  biaxial  loading  on  unstiffened  panels.  If  elastic  anaLysis  is 
considered,  the  influence  of  biaxial  loading  is  known  to  be  slight.  However, 
the  Influence  of  biaxial  loading  may  not  be  negligible  if  plastic  deformation 
is  considered.  In  Reference  4,  the  infLuence  of  biaxial  loading  on  plastic 
stress  intensity  factors  was  examined  (plastic  stress  intensity  factors  are 
related  to  J),  It  was  revealed  that  the  plastic  stress  intensity  factors  are 
reduced  by  appLying  a positive  biaxiaL  load,  i.e,,  the  load  applied  in  two 
directions  is  tensile.  On  the  other  hand,  if  the  biaxial  Load  ratio  is  negative 
(i.e,,  the  Load  parallel  to  the  crack  plane  is  compressive)  the  plastic 
stress  intensity  factors  increase.  The  magnitude  of  the  influence  of  biaxial 
loading  was  shown  to  be  dependent  on  the  strain  hardening  coefficient  of  the 
ma  terial . 

The  influence  of  the  biaxial  loading  on  cracked,  stiffened  panels  was 
studied  in  Reference  3.  This  data  showed  that  the  elastic  stress  intensity 
factors  increased  due  to  the  application  of  positive  biaxiaL  loads  (i.e., 
the  loads  normal  and  perpendicular  to  the  cracked  plane  being  tensile);  and 
that  the  influence  of  biaxial  loading  was  not  negligible.  This  behavior  is 
contrary  to  the  observation  in  the  case  of  unstiffened  sheets.  In  the  present 
study  the  influence  of  biaxial  loading  on  stiffened  sheets  has  been  considered 
using  both  elastic  and  elastic -plastic  analysis.  The  analytical  resuLts  have 
been  compared  with  experimental  data  and  the  influence  of  biaxial  loading  on 
load  transfer  effects  has  also  been  examined. 
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3.2  DESCRIPTION  OF  TEST  PANEL! S) 


Two  identical  test  specimens  were  fabricated  from  the  drawing  of  Figure 

58.  The  skin  material  was  0.063  inch  2024-T3  sheet,  zee  frames  were  0.063- 

inch  thick  2024-T3;  and  the  hat  section  longerons  were  formed  from  7075-T6 
nominally  0.083-inch  thick.  In  the  test  area  the  longerons  were  fastened 
to  the  skin  by  rivets.  It  can  be  noted  from  Figure  58  that  a -1  and  -11 

assembly  are  mentioned.  The  test  panel  identification  is  2-2  and  2-3  to 

reflect  the  skin  material  identification  code.  Between  frames  the  four  bay 
test  area  contained  a centrally  located  initial  slot  in  the  skin  of  4.0-inch 
half  length  and  8.0-inch  overall  length. 

Table  XI  lists  the  materials,  cross  sectional  areas,  and  other  pertinent 
information  for  these  two  panels. 


TABLE  XI  FRAME -LONGERON  STIFFENED  PANEL  AREAS  AND  MATERIAL 

IDENTIFICATION 


PANEL 

NUMBER 

MEASURED  THICKNESS 
( Inches) 

TOTAL  CROSS  SECTIONAL 
AREAS  (Inches^) 

SKIN 

MATERIAL/ 
AND  I.D. 
(See 

Ref.  2) 

PERCENT- 
AGE OF 
TOTAL 
LONGERON 
AREA 

TO  SKIN 

F rame 

Long- 

eron 

Skin 

F rame 

Long- 

eron 

Skin 

2-2 

.063 

.083 

.0635 

.94 

2.3535 

2.540 

2024 -T3/ 
2-2 

92.7 

Total  Panel  Cross 
Sectional  Area  = 4.8935 

2-3 

.063 

.083 

.0635 

.94 

2.3535 

2.540 

2024 -T3/ 
2-3 

92.7 

Total  Panel  Cross 
Sectional  Area  = 4.8935 

A description  of  the  loading  arrangement  will  follow  in  Section 
3.4.  Twenty  strain  gages  were  placed  at  various  locations  on  the  skin, 
longeronls),  and  frame(s).  A standard  beam  type  clip  gage  was  mounted  on  the 
skin  side  of  both  panels  to  record  crack  opening  displacement  at  the  crack 
and  panel  centerline. 
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3,3  FINITE  ELEMENT  MODEL  AND  ANALYSIS 


The  analysis  of  the  biaxial ly  loaded  panel  was  conducted  using  two 
different  biaxial  load  ratios*  The  finite  element  modeling*  elastic*  and 
elastic -plastic  analyses  of  the  panel  are  discussed  in  the  following 
subsections. 

3.3. L Finite  Element  Modeling  of  Biaxial  I y Loaded  Panels 

The  finite  element  modeling  used  is  similar  to  that  used  in  the  analysis 
of  the  Phase  II  panels  (Reference  1).  The  panel  is  essentially  modeled  as  a 
two-dimensional  structure*  Triangular  membrane  elements  were  used  to  model 
the  skin . The  longerons  of  this  panel  are  modeled  like  the  stiffeners  of  the 
wing  channel  panel  of  Phase  II  (Section  3*1  of  Reference  i>*  The  connected 
leg  of  the  longeron  is  modeled  as  a rod  element  and  the  two  upstanding 
portions  of  the  longeron  (web  portions)  are  combined  to  represent  rectangular 
membrane  elements.  The  thickness  of  the  membrane  elements  is  equal  to  the 
combined  thickness  of  the  two  outstanding  portions  (web  elements).  The  two 
flange  portions  and  lips  are  combined  to  represent  one  rod  element.  The 
frame  is  modeled  in  a similar  way  as  the  longeron.  The  connected  portion  is 
modeled  as  a rod  element  and  the  outstanding  leg  as  a membrane  element.  The 
flange  and  lip  portion  of  the  frame  were  combined  to  represent  one  rod 
element  * 

The  portion  of  the  shear  clip  connected  to  the  skin  was  modeled  as  a rod 
element  and  that  portion  connected  to  the  frame  was  modeled  as  a membrane 
element.  The  rivets  connecting  the  skin  to  the  shear  clips  were  modeled 
using  the  flexible  fastener  model.  The  first  nine  rivets  from  the  crack  plane, 
connecting  the  skin  to  the  longeron,  were  modeled  as  individual  shear  elements* 
These  shear  elements  were  proportioned  to  represent  a flexible  fastener  model. 
The  portion  of  the  panel  beyond  the  first  nine  rivets  was  assumed  to  be 
connected  by  continuous  shear  elements  as  for  the  angle  stiffened  panels  (refer 
to  Subsection  2. 2.2.1).  The  shear  elements  were  proportioned  to  account  for 
fastener  flexibility.  Only  one  quarter  of  the  panel  was  modeLed  for  finite 
element  analysis.  The  finite  element  model  of  the  biaxial  panel  is  shown  in 
Figure  59. 

3.3,2  Elastic  and  Elastic -F lastic  Analysis  of  the  Biaxial  Panel 

The  elastic  and  elastic -plastic  analysis  of  the  biaxial ly  loaded  panel 
was  performed  to  study  the  influence  of  biaxial  loading  on  load  transfer  to 
stringers,  crack  openings,  and  J values.  The  elastic  and  elastic -plastic 
analysis  of  this  panel  is  discussed  in  the  following  subsections, 

3,3.2. 1 Elastic  Analysis  of  Biaxial  Panel 

The  elastic  analysis  of  the  biaxial  paneL  was  performed  for  various  crack 
lengths  (a  = 4.5,  5*75,  8.00,  9.00,  11,0  inches)  with  a uniaxial  to  biaxial 
ratio  of  1:0  (i,e.,  no  biaxial  load  applied).  Two  different  contours  were 
used  to  evaluate  /J  for  a half  crack  length  of  4.5  inches.  These  contours 
are  shown  In  Figure  60.  The  /J  values  determined  for  these  two  contours 
were  within  one  percent  of  each  other.  Thus  for  other  c rack  lengths  only 
Contour  I (Figure  60)  was  used  to  evaluate  the  square  root  of  J values. 
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Figure  59 . Finite  Element  Model  of  Biaxial  PaneL(s) 


o 
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Figure  60.  Contours  Used  to  Evaluate  J for  Biaxial  Panels 


Values  of  square  root  J obtained  for  various  crack  lengths  with  the  1:0 
biaxial  ratio  are  shown  in  Figure  61.  The  /T  value  increases  with  crack  Length. 
When  the  crack  tip  is  in  the  vicinity  of  the  next  stringer  (longeron) 
the  values  decrease.  Once  the  crack  tip  is  beyond  the  centerline  of  the 
Longeron,  the  square  root  J values  increase  again.  The  stresses  in  the 
central  stringer  (longeron)  for  various  crack  lengths  are  indicated  in 
Figure  62.  For  the  same  appLled  stress,  the  stress  in  the  central  stringer 
increases  with  crack  length. 

An  elastic  analysis  of  the  biaxial  panel  was  also  performed  at  a biaxiaL 
load  ratio  of  L:i  for  a half  crack  length  of  4.5  inches.  For  this  c rack 
Length,  the  plot  of  /J  versus  applied  stress  is  shown  in  Figure  63,  for  a 
biaxial  ratio  of  1:L  and  also  1:0  (no  biaxial  load  applied).  The  biaxial 
load  ratio  1:1  here  refers  to  ratio  of  total  load  applied  in  y direction 
to  load  applied  in  x direction,  Lt  is  surprising  to  note  that  for  the 
same  stress  applied  in  the  uniaxial  or  y direction  (p)  the  square  root  J 
value  increases  due  to  the  application  of  a biaxiaL  load.  The  values  due  to 
the  application  of  a biaxiaL  load  being  L3  percent  higher  than  those  due  to  a 
uniaxial  load  only.  This  behavior  is  contrary  to  that  observed  in  unstiffened 
panels  where  it  is  known  that  the  effect  of  biaxial  tensile  load  with  the 
same  sign  as  the  uniaxial  Load,  is'  to  reduce  the  crack  opening  displacement 
and  stress  intensity  factors  (Reference  4),  For  unstiffened  panels  the  effect 
of  biaxiaL  load  on  elastic  analysis  is  known  to  be  of  a second  order  of 
magnitude  (Reference  5),  As  J is  related  to  the  stress  intensity  factor, 
this  increase  of  square  root  of  J with  application  of  biaxial  load  for 
stiffened  panels  indicates  that  the  stress  intensity  factor  ahead  of  the  crack 
tip  increases  due  to  biaxial  load.  Similar  increases  in  stress  intensity 
factors  due  to  biaxial  Loads  for  stiffened  panels  were  observed  in  Reference  3. 
This  is  due  to  the  fact  that  biaxial  load  causes  less  load  transfer  to  the 
central  stringer  and  consequently  larger  crack  surface  openings  or  displacements. 
The  crack  surface  openings  are  related  to  the  stress  intensity  factors  and 
hence, the  stress  intensity  factors  increase  due  to  the  application  of  a biaxial 
load  in  a stiffened  panel.  Figure  64  shows  the  plot  of  stress  in  the  central 
stringer  (longeron),  at  the  plane  of  the  crack,  for  biaxial  load  ratios  of  1:1 
and  1:0,  The  application  of  a 1:1  biaxial  load  causes  a 13  percent  reduction 
in  stringer  or  Longeron  stress.  The  crack  surface  openings  as  a function  of 
distance  from  the  center  of  the  crack  (along  the  crack  Line)  , for  uniaxial 
load  (biaxial  load  equals  zero)  and  biaxial  load  ratio  1:1  are  shown  in 
Figure  65  • The  crack  surface  openings  or  displacements  at  the  center  of  crack 
increase  L4  percent  over  the  uniaxial  load  case  (biaxial  ratio  1:0)  when  a 1:1 
biaxial  load  ratio  is  applied.  Thus  for  a stiffened  paneL  the  elastic  analysis 
indicates  that  the  effect  of  biaxial  load  is  to  increase  crack  opening  and  J 
and  reduce  load  transfer  to  the  central  stringer.  The  influence  for  this  panel 
geometry  was  determined  to  be  between  13  to  15  percent , 

The  stresses  in  the  stringers  ( longerons)  located  at  9 inches  and  18 
inches  from  the  centerline  of  the  panel  and  crack  are  shown  in  Figure  66,  The 
stresses  in  the  stringer,  9 inches  from  the  centerline  of  the  crack  are  reduced 
by  23  percent  and  the  stresses  in  the  stringer  18  inches  away  are  reduced  by 
4 percent  due  to  a biaxial  load  ratio  of  1:1. 
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Figure  61.  Square  Root  of  J for  Varying  Crack  Lengths  as  a Function  of  Applied  Stress 
Biaxial  Panel 
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Figure  62.  Stress  in  the  Central  Longeron  as  a Function  of  Gross  Area 
Panel  Stress  for  Varying  Crack  Length  - Biaxial  Panel 
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Square  Root  of  J , /J 


Figure  63.  Elastic  /j  Values  as  a Function  of  Applied  Stress  for  Two  Biaxial  Load  Ratios 
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Figure  64,  Stress  in  the  Central  Longeron  as  a Function  of  Applied  Stress  for 
Two  Biaxial  Load  Ratios 
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Distance  From  Center  of  the  Crack  and  Panel,  (Inches) 

Figure  65.  Total  Crack  Surface  Displacements  as  a Function  of  Distance  Along  the 
Crack  for  Two  Biaxial  Load  Ratios 
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Stress  in  Longerons,  (ksi) 


v 


Figure  66.  Stresses  In  Adjacent  Longerons  for  Two  Biaxial  Load  Ratios 
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3. 3. 2. 2 Elastic-Plastic  Analysis  of  Biaxial  Panel 


The  elastic-plastic  analysis  of  the  biaxially  loaded  panel  based  on 
Dugdale  type  strip  plastic  zone  could  not  be  performed  under  biaxial  loading 
due  to  the  fact  that  the  Dugdale  model  is  only  applicable  to  uniaxial  load 
cases.  The  application  of  a biaxial  stress  will  change  the  state  of  the 
stress  and  hence  the  plastic  zone  ahead  of  the  crack  tip.  The  size  of  the 
Dugdale  plastic  zone  is  determined  by  removing  stress  singularities  ahead  of 
the  crack.  The  effect  of  biaxial  loading  on  stress  singularities  is  of 
second  order  and  hence  the  influence  of  biaxial  loading  on  a Dugdale  type 
plastic  zone  will  be  very  small.  Thus  a Dugdale  analysis  in  biaxial  loading 
will  not  give  a realistic  representation  of  plasticity  ahead  of  the  crack  tip. 
Considering  these  factors  it  was  felt  more  appropriate  to  perform  the  elastic- 
plastic  analysis  of  the  biaxial  panel  based  on  P rand t 1 - Reuss  material  behavior. 
The  elastic- plastic  analysis  based  on  P rand 1 1 - Reuss  material  behavior  for 
various  load  increments  involves  considerable  computer  time.  Thus,  it  was 
decided  to  perform  an  elastic- plastic  analysis  assuming  Prandt 1 - Reuss  material 
behavior  for  only  one  crack  length.  The  length  selected  was  8.30  inches 
or  the  crack  length  of  the  tested  panels  with  biaxial  load  ratios  of  1:0, 

1:1/3,  1:2/3,  and  1:1.  The  corresponding  applied  load  increments  were  10,  20, 
30,  40,  and  45  ksi.  The  J integral  values  were  computed  using  Contour  1 
shown  in  Figure  60  . The  variation  of  /J  with  applied  stress  for  a half  crack 
length  of  4.25  inches  and  various  biaxial  ratios  is  shown  in  Figure  67.  For 
the  same  applied  load  normal  to  the  crack  (p)  the  /J  values  increase  as  the 
load  in  the  x direction  is  increased,  i.e.,  biaxial  load  ratio  is  increased. 
This  behavior  is  similar  to  that  observed  for  the  elastic  analysis.  The  stress 
at  which  nonlinear  effects  are  observed  decreases  as  the  biaxial  load  ratio 
increases.  The  variation  of  stress  in  the  central  stringer,  at  the  plane  of 
the  crack,  as  a function  of  applied  stress,  for  various  biaxial  load  ratios 
is  shown  in  Figure  68  . Here  again  the  behavior  is  similar  to  the  elastic 
analysis  where  the  stress  in  the  central  stringer  decreases  as  the  biaxial 
load  ratio  increases,  i.e.,  load  in  the  x direction  is  increased. 

The  influence  of  biaxial  loading  on  the  plastic  zone  size  ahead  of  the 
crack  tip  is  shown  in  Figure  69  . It  is  seen  that  at  an  applied  stress 
(normal  to  the  crack)  of  40  ksi,  the  plastic  zone  size  decreases  with  the 
increase  in  the  biaxial  load  (the  applied  load  in  x direction  being  of  same 
sign  as  the  load  normal  to  the  crack).  Similar  reduction  in  the  size  of  the 
plastic  zone  due  to  application  of  biaxial  load  was  also  observed  in  the 
analysis  of  Reference  4. 

3.4  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  details  of  load  application,  strain  gage  location,  strain  surveys, 
etc.,  were  common  to  both  the  2-2  and  2-3  panels  and  a general  description 
of  each  will  follow.  All  instrumentation  employed  to  record  and  analyze 
data,  i.e.,  load,  COD,  strain  were  identical  to  that  used  for  the  uniaxiaily 
loaded  wing  panels.  A description  of  this  test  and  data  reduction  arrange- 
ment is  given  in  Section  2.2.3. 1 and  will  not  be  repeated  here. 
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Figure  67.  Square  Root  of  J for  Varying  Biaxial  Load  Ratios  Using  Prandtl-Reuss 
Material  Assumptions 
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Figure  68.  Stress  tn  Central  Longeron  for  Varying  BiaxiaL  Load  Ratios  Using 
Prandti-Reuss  Material  Assumptions 
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Figure  69.  Plastic  Zone  Sizes  for  Three  Biaxial  Load  Ratios  Using  Prandtl-Reuss 
Material  Assumptions 
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The  lateral  loading  was  applied  through  two  12),  ten- inch  hydraulic 
cylinders  each  capable  of  providing  a 120,000  pound  load  through  a whiffle- 
tree  arrangement  connected  to  the  frames/skin  (see  Section  B-B  of  Figure  58) 
on  both  sides  of  the  panel*  A view  of  the  overall  structural  test  setup 
is  shown  in  Figures  70  and  71  in  which  the  support  frame  for  the  lateral 
load  cylinders  and  whiffle tree  can  be  seen  in  the  500,000  pound  load  frame. 

A four  channel  EDISON  load  maintainer  was  used  to  control  load  for  these  side 
cylinders. 

It  will  be  noted  that  in  the  loading  arrangement  shown  in  Figure  70 
that  some  panel  bending  is  unavoidable  since  the  lateral  l side)  loading  is 
fixed  and  not  free  to  move  in  a vertical  direction  when  load  is  applied  along 
the  panel  length.  As  with  the  previous  wing  panels  care  was  taken  to  minimize 
inplane  bending  by  padding  at  the  grip  ends.  It  will  be  noted  in  the  dis- 
cussions of  the  test  data  for  the  individual  panels  that  the  lateral  con- 
straint presented  a rotational  bending  problem  ooly  at  the  highest  lateral 
side  loads, 

3,4,1  Fatigue  PrecrackinE,  Strain  Data  and  Fracture  of  Biaxial  Loaded  Paneljj 

To  start  from  a natural  crack  condition  in  the  two  panels  each  was  fatigue 
precracked  approximately  one  quarter  inch  from  the  ends  of  the  machined  slot 
using  tension- tension,  sinusoidal  loading.  An  approximation  of  the  stress 
intensity  for  this  panel  was  made  using  the  solutions  of  Reference  3.  Assuming 
a maximum,  cyclic  stress  intensity  of  12  ksi  J inch  in  both  cases.  It  was 
subsequently  determined  that  the  stress  intensity  was  conservative  and  a 
change  was  made  to  higher  load  levels.  The  fatigue  stress  ratio  was  maintained 
at  between  0.1  and  0.25  to  restrain  local  crack  bending  at  lower  load  levels. 
Table  XII  summarizes  the  fatigue  loadings  and  total  crack  lengths  for  the  two 
biaxial  panels.  Measurement  of  final  fatigue  crack  length  was  accomplished 
after  panel  failure. 


TABLE  XU  FATIGUE  PRECRACKING  AND  CRACK  LENGTH  DATA,  2-2  AND 
2-3  BIAXIAL  PANELS 


PANEL 

FATIGUE  LOADS 

NUMBER  OF 

CYCLES 

FATIGUE  CRACK  LENGTH 
( Inches) 

MAXIMUM 

Ikips) 

MINIMUM 

ikips) 

LEFT 

RIGHT 

TOTAL 

2-2 

19.5 

2.0 

31,000 

- 

- 

- 

24.0 

2.4 

44,000 

*06 

.09 

8. 15 

28,0 

7.0 

11,800 

.20 

.20 

8.40 

2-3 

28  *0 

7.0 

22,000 

.26 

.24 

8.50 
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Figure  70.  Stiffener-Side  View  of  Loading  Arrangement  for  Biaxial  Panels 
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Figure  71.  Side  View  of  Lateral  Loading  Arrangement  for  Biaxial  Panels 
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Figures  72  and  73  show  the  front  (skin)  side  view  of  the  biaxial  panel 
in  the  test  machine.  It  will  be  noted  that  the  diamond  shaped  doublers  are 
used  to  aid  in  transmitting  load  to  the  skin  from  the  whiff elt reel s) . The' 
test  section  is  that  area  shown  in  Figure  73  between  frames.  Also  visible  is 
the  starter  slot  in  the  skin  in  Figure  72.  The  positioning  of  the  displacement 
gage  points  and  most  strain  gage  locations  can  be  noted  on  both  views. 

After  fatigue  cycling  and  zero  balancing  of  all  strain  gages  and  dis- 
placement gage  with  the  panel  in  a free  hanging  position  a series  of  strain/ 
load  survey  was  accomplished  on  each  panel  by  loading  to  specific  values  of 
biaxial  ratio.  The  biaxial  ratio  here  is  defined  as  ratio  of  lateral  (frame)  or 
side  load  to  axial  (lengthwise)  load.  Care  was  taken  to  perform  this  survey 
within  the  elastic  range  for  each  panel.  The  results  of  these  surveys  and 
strains  during  the  fracture  and  slow  tear  portion  of  the  loading  cycle  are 
given  in  Table  XLL1  for  the  2-2  panel  and  Table  XIV  for  the  2-3  panel.  The 
reference  location  of  each  strain  gage  is  shown  in  Figure  1U  located  between 
Tables  XIII  and  XIV. 

Load  versus  total  displacement  gage  output  are  shown  in  Figures  75  and 
76  for  the  2-2  panel  and  in  Figure  77  for  the  2-3  panel.  It  is  of  interest 
to  note  the  increase  in  crack  opening  displacement  with  increasing  lateral 
load  under  constant  axial  load  (see  Figure  75)*  This  behavior  was  predicted 
by  the  analytical  results.  The  loading  to  fracture  of  the  2-2  panel  was  of 
ramp  loading  (100,000  lbs. /min.)  to  60  kip  lateral  load,  keeping  that  load 
constant  and  increasing  axial  load.  However,  loading  was  interrupted  due  to 
an  equipment  malfunction  during  the  fracture  loading  process.  This  resulted 
in  a stable  tear  to  an  overall  crack  length  of  11.28  inches  (5.67  inches  on  the 
left  side  and  5.61  inches  on  the  right  side)  at  which  time  the  load  was  removed 
and  a subsequent  loading  exercise  was  attempted  (see  Figure  75).  At  this 
crack  length  the  panel  was  once  again  ramp  loaded  at  100,000  pounds  per  minute 
in  both  axial  and  lateral  load.  The  lateral  loading  was  stopped  at  60  kips 
and  the  axial  load  increased  at  the  same  rate.  Failure  occurred  in  this  2-2 
panel  at  195.1  kips. 

After  the  initial  strain  survey  in  the  2-3  panel  (see  Table  XIV)  the 
panel  was  loaded  at  100,000  pounds  per  minute  to  a lateral  and  axial  load  of 
120,000  pounds.  At  that  point  axial  loading  was  increased  to  failure  at 
181  kips. 

Motion  picture  coverage  was  initiated  during  the  slow  tear  and  fracture 
process  for  both  biaxially  loaded  panels.  In  both  cases  the  tear  to  fracture 
and  fracture  were  recorded.  The  failure  sequence  for  each  panel  is  shown  in 
Figures  78  and  79.  It  is  obvious  that  the  2-2  panel  (Figure  78)  failed  in 
the  test  section  and  through  the  longeron/skin  area.  However,  examination  of 
the  failure  sequence  in  Figure  79  for  the  2-3  panel  indicates  that  fracture 
initiated  in  the  area  of  the  attached  doubler  near  the  frame  section  (see 
Frame  No.  9).  The  resulting  failure  of  this  panel  was  explosive.  The  design 
of  the  frame/skin  lateral  loading  attachments,  etc.,  was  placed  at  approxi- 
mately 130  kips  (with  no  margin  of  safety)  therefore  a failure  in  this  area 
could  have  been  avoided  by  employing  a lesser  lateral  load.  However,  in  this 
phase  of  the  program,  with  two  identical  specimens  an  order  of  difference  in 
lateral  loading  between  the  two  specimens  was  the  logical  scheduled  test  plan. 
Therefore,  a constant  60  kip  and  120  kip  side  load  were  selected  as  the  fracture 
conditions  prior  to  test. 
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Figure  72.  Skln-Slde  View  of  Biaxial  Test  Panel 
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Figure  73.  Longeron  and  Frame  Side  View  of  Biaxial  Test  Panel 
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TABLE  XLlI  STRAIN  GAGE  DATA,  2-2  BIAXIAL  PANEL 
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Figure  74.  Strain  Gage  Locations  - Biaxial  Panels 


TABLE  XIV  STRAIN  GAGE  DATA,  2-3  BIAXIAL  PANEL 
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Figure  75.  Load -Di  splacement  Curves  for  Biaxial  Panel  2-2  During  Strain  Survey  and  Slow  Tear 
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Figure  78.  Fracture  Sequence  for  2-2  Biaxial ly  Loaded  Fanel  (60  kips  Lateral 
Load ) 
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Figure  78.  Fracture  Sequence  for  2-2  Biaxially  Loaded  Panel  (60  kips  Lateral 
Load)  (Continued) 
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Figure  79.  Fracture  Sequence  for  2-3  Biaxially  Loaded  Panel  (120  kips  Lateral 
Load ) 
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Figure  79.  Fracture  Sequence  for  2-3  Biaxial ly  Loaded  Panel  (120  kips  Lateral 
Load)  (Continued) 
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Fracture  surfaces  of  the  2-2  panel.  (60  kip  lateral  load  at  fracture)  are 
shown  in  Figure  80  for  the  skin  side  and  in  Figure  81  for  the  longe ron/ f rame 
side.  At  fracture  the  crack  ran  to  (and  through)  the  fastener  holes  in  the 
adjacent  longerons.  The  frames  remained  intact  and  four  of  the  longerons 
fractured  along  the  crack  path.  In  the  longeron  which  did  not  fail  in  the 
test  section  shearing  of  rivets  was  quite  prevalent  (see  Figure  81). 

Figures  82  and  83  show  the  fracture  surfaces  of  the  skin  and  angle  side 
for  the  2-3  panel  (120  kip  lateral  load  at  failure).  As  indicated  in  the 
failure  sequence  for  the  2-3  panel  in  Figure  79  the  path  was  jagged  and 
failure  progressed  from  one  doubler  across  the  panel  to  the  unzippering  in 
the  longeron  area  on  the  opposite  side  of  the  panel. 

3. A. 2 Comparison  of  Experimental  and  Analytical  Data 

Figure  74  shows  the  location  of  strain  gages  for  the  biaxial  panels. 
Tabulated  in  Tables  XIII  and  XIV  are  the  strain  gage  readings  for  various 
applied  loads.  Figure  84  shows  the  variation  of  strain  in  the  central 
stringer  with  applied  stress  for  various  biaxial  stress  ratios.  The  strains 
in  the  stringer  at  1.3  and  2.5  inches  from  the  crack  plane  show  good  correla- 
tion with  analytical  results  for  biaxial  ratios  of  1:0,  1:0.3  and  1:1. 

However  the  strain  at  0.6  inch  from  the  crack  plane  differs  considerably  from 
the  analytical  results.  Figure  85  shows  the  variation  of  strain  in  the 
central  stringer  (longeron)  and  the  stringer  9 inches  away  from  the  centerline 
of  the  crack,  with  applied  stress  for  a constant  lateral  load  (Px)  of  60  kips. 

This  constant  lateral  load  of  60  kips  results  in  different  biaxial  load  ratios, 
depending  on  the  applied  load  normal  to  the  crack  surface.  The  experimental 
strains  are  shown  for  an  initial  crack  length  of  8.4  inches  and  "post  tear" 
crack  length  of  11.28  inches.  Analytical  strains  for  a crack  length  of  8.5 
inches  are  also  noted  in  Figure  85.  The  correlation  between  experimental  and 
analytical  strains  in  the  stringer  (longeron)  located  9 inches  from  the 
centerline  of  the  crack,  is  good.  The  strains  in  the  central  stringer  show  a 
good  agreement  with  analytical  results  at  all  but  low  applied  stresses.  The 
analytical  strains  for  a crack  length  of  11.28  inches  were  not  calculated  so 
no  comparison  could  be  made. 

Figure  86  shows  the  variation  of  strains  in  the  frame  members  for  three 
biaxial  stress  ratios,  1:0,  1:0.5  and  1:1.  Analytical  strains  obtained  from 
the  finite  element  analysis  are  also  noted.  For  zero  biaxial  load  the  strains 
in  the  flange  and  upstanding  leg  (web)  show  a good  agreement  with  analytical 
results.  However,  strains  in  the  attached  leg  show  a considerable  variation 
from  the  analytical  results.  At  other  biaxial  load  ratios,  the  strains  in 
the  flange  portion  of  the  frame  show  good  agreement  with  the  analytical  data 
but  strains  in  the  attached  and  outstanding  legs  (web)  differ  considerably 
from  the  analytical  results.  This  discrepancy  between  analytical  and  experi- 
mental attached  leg  strains  is  perhaps  due  to  the  simplified  attachment  model 
combined  with  the  method  of  load  application  used  in  the  finite  element 
analysis.  In  the  finite  element  model  individual  rivets  connecting  the  frame 
to  the  skin  were  not  modeled  but  were  represented  by  shear  elements  having  a 
width  of  one  inch.  It  was  also  assumed  that  at  the  side  of  the  panel  the 
biaxial  (lateral)  load  was  uniformly  distributed  over  a width  of  13  inches. 
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Figure  80.  View  of  Fracture  Surface  from  Skin  Side,  2-2  Biaxial  ly  Loaded  Panel 
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Figure  81.  View  of  Fracture  Surface  from  Longeron/F  rame  Side,  2-2  Biaxialiy  Loaded  Panel 
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Figure  82.  View  of  Fracture  Surface  from  Skin  Side,  2-3  Biaxially  Loaded  Panel 
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Figure  83.  View  of  Fracture  Surface  from  Longeron /F rame  Side,  2-3  Biaxially  Loaded  Panel 
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Figure  84.  Strain  in  CentraL  Longeron  Away  From  the  Crack  at  Three  Biaxial 
Load  Ratios 
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Figure  85.  Strain  in  Longeron(s)  at  0.6  Inches  Away  from  the  Crack, 
Biaxial  Panel  2-2 
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Strain  In  Frame,  (Inches/Inch  x 10^) 


Figure  86.  Strain  In  Frame  Member  as  a Function  of  Applied  Stress  for 
Three  Biaxial  Load  Ratios 
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Figure  87  shows  the  experimentally  observed  crack  surface  openings  at 
the  center  of  crack  (and  panel)  as  a function  of  applied  stress  for  zero 
biaxial  load  and  a biaxial  load  ratio  of  1:1*  Also  indicated  In  Figure  H7 
are  the  analytically  derived  crack  openings.  The  agreement  between  analy- 
tical and  experimental  data  is  excellent.  The  application  of  a biaxial 
(lateral)  load  increases  the  crack  openings  as  was  observed  both  analytically 
and  experimentally  for  these  panels. 

3.5  RESIDUAL  STRENGTH  PREDICTION 


For  uniaxial ly  loaded  panels  the  residual  strength  prediction  was  based 
on  a DugdaLe  type  elastic -plastic  analysis  and  /Jr  resistance  curve  of  the 
material*  The  /Jp  resistance  curve  of  the  material  was  obtained  from  tests 

on  the  crack  line  wedge  loaded  specimen.  There  are  no  biaxial  loads  applied 
in  obtaining  /Jr  resistance  curve.  The  influence  of  biaxial  loading  on  the 

resistance  curves  of  a material  has  not  been  studied*  It  was  not  feasible  to 
study  this  influence  in  this  program,  hence  it  was  decided  to  use  the  /J^ 

resistance  curve  of  the  material  without  biaxial  effects*  As  discussed 
earlier,  Dugdale  type  elastic-plastic  analysis  could  not  be  performed  for 
biaxial ly  loaded  panels  and  the  realistic  analysis  of  this  panej.  should  be 
based  on  F rand t l - Reuss  material  behavior*  In  order  to  obtain  / J versus  crack 
Length  curves  for  various  applied  loads,  the  elastic-plastic  analysis,  assuming 
P randt l - Reuss  material  behavior,  would  have  to  be  performed  for  various  crack 
lengths  and  applied  loads.  This  would  require  extremely  large  computer  run 
times  with  prohibitive  cost.  Ln  the  present  study  an  e las t ic - plasti c analysis 
was  performed  for  only  one  crack'  Length  and  elastic  analysis  was  performed  for 
various  crack  lengths*  In  order  to  obtain  /J  versus  a curves,  certain  approxi- 
mations were  made.  It  was  assumed  that  for  a fixed  applied  stress,  the  per- 
centage increase  in  /J  due  to  biaxial  load  and  P randt l - Reuss  material  behavior 
over  elastic  values  was  the  same  for  all  crack  lengths.  From  Figure  67,  the 
/J  value  for  1:1/3  biaxial  load  is  10  percent  higher  than  that  for  the  elastic 
1:0  load,  at  an  appiied  stress  of  40  ksi.  Thus  it  is  assumed  that  the  /J 
values  for  1:1/3  biaxial  load  will  be  IQ  percent  higher  than  those  /I  values 
based  on  an  elastic  analysis  for  all  crack  lengths  at  an  applied  stress  of 
40  ksi.  The  percentage  increase  in  /J  value  for  other  applied  stresses  will  be 
different.  With  this  assumption  / J versus  a curves  for  various  applied  stresses 
are  shown  in  Figure  88  for  a biaxial  Load  ratio  of  1:1/3* 

The  stresses  in  the  central  stringer  of  the  biaxial  panel  for  various 
applied  stresses  based  on  P rand t i - Reuss  material  behavior  are  shown  in 
Figure  68*  The  uLtimate  stress  for  the  stringer  material  from  tensile  test 
data  was  found  to  be  80  ksi.  Using  this  ultimate  stress  and  the  stringer 
stress  for  a biaxial  load  ratio  of  1:1/3  the  failure  stress  is  predicted  at  an 
applied  stress  of  56,5  ksi*  In  Figure  88  the  resistance  curve  of  the  skin 
materiaL  l 2024-T3, 0,063"  (LT))  has  been  pLotted  at  a half  crack  length  of  4,25 
inches.  There  is  considerable  slow  tear  at  an  appLied  stress  of  40  ksi.  Thus 
the  stringer  stresses  shown  in  Figure  68  will  not  be  valid*  Also  from 
Figure  88  it  is  seen  that  there  is  no  point  of  instability  within  the  portion 
of  resistance  curve  of  the  material  obtained  from  the  CLWL  tests. 
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Figure  87.  Total  Crack  Opening  Displacements  for  Two  Biaxial  Load  Ratios 
as  a Function  of  Applied  Stress 
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Figure  88,  Square  Root  o£  J Versus  Crack  Length  Curve  for  Biaxial  Panel 
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However,  if  the  resistance  curve  of  the  material  is  extrapolated  the  point 
of  instability  occurs  at  an  applied  stress  of  42  ksi.  The  2-2  panel  slowly 
tore  to  a half  crack  length  of  approximately  5.63  inches  and  then  it  was 
unloaded.  The  2-2  panel  was  tested  again  after  a delay  period  due  to  equipment 
malfunction  with  this  same  crack  length.  The  failure  load  of  the  panel  was  195.1 
kips  (see  Table  XIII)  and  there  was  no  crack  arrest  in  the  panel.  In  Figure  88, 
replotting  the  resistance  curve  C/Jr^)  of  the  material  at  a half  crack  length  of 

5.63  inches  shows  that  the  point  of  instability  occurs  at  an  applied  stress  of 
40  ksi  and  /J  values  of  the  panel  are  higher  than  the  resistance  curve  of  the 
material.  Hence,  there  should  be  no  crack  arrest  in  this  2-2  panel  at  this 
crack  length.  The  predicted  failure  load  corresponding  to  a stress  of  40  ksi 
is  195.740.  This  is  within  less  than  one  percent  of  the  actual  failure  load. 

It  must  be  emphasized  that  the  /J  versus  a curves  in  Figure  88  are  only 
approximate  and  the  /Jr  curve  is  based  on  uniaxial  data  and  hence  the  pre- 
dicted load  is  only  an  approximation. 

The  approximated  square  root  of  /J  versus  a curves  for  the  biaxial  load 
of  1:2/3  were  not  plotted.  Referring  to  Figure  67,  it  can  be  seen  that  for 
a 1:2/3  biaxial  load  ratio  the  /J  values  differ  considerably  from  the  elastic 
values  and  hence  the  assumptions  made  in  obtaining  the  /J  versus  a curves 
in  biaxial  load  ratio  1:1/3  cannot  be  justified.  In  the  2-3  panel,  tested 
at  a biaxial  load  ratio  of  1:2/3,  the  crack  tore  to  the  rivet  hole  at  the 
adjoining  stringer  and  was  arrested.  With  subsequent  increase  in  load,  the 
central  stringer  fractured  and  the  failure  started  at  the  doublers  through 
which  biaxial  load  was  applied  to  the  skin.  Figure  62  shows  the  stresses  in  the 
central  stringer  based  on  elastic  analysis  and  zero  biaxial  load.  Using  the 
measured  ultimate  stress  of  the  stringer  material  (80  ksi)  and  the  stress  in  the 
stringer  for  a half  crack  length  of  9 inches  (crack  length  at  which  the  crack 
was  arrested  at  the  longeron  rivet  hole)  the  failure  of  the  stringer  is 
predicted  as  40.7  ksi.  Due  to  biaxial  loading  the  stresses  in  the  stringer 
will  be  smaller  than  those  shown  for  zero  biaxial  load.  If  an  elastic- 
plastic  analysis  is  considered  the  stresses  will  be  larger  than  those  shown 
in  Figure  62.  Thus,  the  stresses  given  in  Figure  62  may  be  taken  as  a 
rough  approximation  of  the  stresses  in  the  stringer  for  the  2-3  panel.  The 
predicted  "approximate"  failure  load  is  40.7  ksi  and  the  actual  failure  load 
of  the  2-3  panel  was  37  ksi. 

In  order  to  predict  the  failure  load  more  accurately  an  e last ic - plas t ic 
analysis  should  be  performed  assuming  P rand 1 1 - Reuss  material  behavior  with 
biaxial  load  for  various  crack  lengths.  From  this  analysis  /J  versus  crack 
length  curves  can  be  plotted  for  various  applied  loads.  These  curves  should 
be  used  along  with  the  /J^  resistance  curve,  obtained  under  biaxial  loading, 

to  predict  the  residual  strength  for  skin  critical  cases. 
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IV  SUMMARY  OF  OVERALL  PROGRAM  AND  PREDICTION  METHOD 


In  an  attempt  to  summarize  the  activity  conducted  during  the  time  span 
of  this  program  it  is  believed  that  by  listing  the  major  tasks,  and  including 
a summary  of  each,  a fuller  understanding  of  the  scope  can  be  appreciated. 
Therefore  each  major  task  will  be  listed  in  this  Section  along  with  the 
significant  findings  of  the  various  phases. 

4.1  PROGRAM  SUMMARY 

The  objective  of  Phase  I of  this  study  was  to  analyze  and  evaluate  the 
state-of-the-art  of  plane  stress  fracture  analysis  and  to  define  the  strong 
and  weak  points  of  those  available  methods. 

It  was  determined  that  the  methods  could  not  adequately  treat  the 
problems  of  slow  stable  tear  (prior  to  fracture)  and  associated  crack  tip 
plasticity,  which  is  prevalent  in  higher  toughness  materials.  The  conclusion 
was  that  an  ideal  residual  strength  prediction  technique  should  account  for 
both  of  these  behaviors  in  the  method  to  be  developed.  With  the  ideal  method 
in  mind  the  development  of  the  analytical  technique  and  proposed  failure 
criterion  were  part  of  the  tasks  of  Phase  II. 

Phase  II  represented  one  level  of  complexity  for  the  developed  residual 
strength  prediction  method  which  was  intended  to  treat  problems  of  cracks  in 
complex,  uniaxially  loaded  structural  components.  To  achieve  this  objective 
crack  growth  resistance  data  were  developed  for  a wide  range  of  aluminum  and 
titanium  alloys.  These  data  were  eventually  used  to  establish  a skin  critical 
fracture  stress  in  the  six  zee  stiffened  panels  examined  during  Phase  11. 

The  six  panels  were  grouped  in  three  sets  of  two  each  containing  the  same 
crack  size  but  three  different  methods  of  attachment  (riveted,  bolted,  and 
bonded).  This  allowed  for  refinement  of  the  fastener  model  and  procedural 
steps.  The  conditions  for  a skin  or  stiffener  critical  situation  were  also 
established.  During  the  development  of  the  analytical  method  the  path  in- 
dependency of  the  J integral  for  stiffened  panels  was  demonstrated  along 
with  the  first  complete  elastic -plastic  analysis  of  a stiffened  structure 
assuming  both  Dugdale  and  P rand t 1 - Reuss  material  behavior.  In  order  to 
establish  an  e las tic  - plastic  data  base,  a nonlinear  analysis  was  performed 
on  the  crack  line  wedge  loaded  specimen  geometry  which  permitted  the  intro- 
duction of  the  Jr  resistance  curve  for  skin  critical  structure.  The  results 

of  several  fracture  criterion  sub-studies  indicated  that  crack  extension 
occurs  at  very  low  values  of  stress  intensity  for  anticipated  plane  stress 
material  behavior  and  the  specimen  independency  of  the  crack  growth  resistance 
curve. 

A ten  step  residual  strength  predictive  procedure  was  established  based 
on  the  results  of  Phase  II  of  the  program. 

The  Phase  III  objective,  reported  on  in  this  report,  was  to  test,  evaluate, 
and  verify  the  capability  of  the  residual  strength  predictive  technique  of 
Phase  11.  This  involved  one  additional  level  of  panel  complexity  and  the 
addition  of  a biaxial  loading  condition  for  fuselage  structure.  During  the 
course  of  the  analysis  of  these  biaxially  loaded  panels  the  effect  of  load 
transfer  to  intact  stringers  was  determined  and  the  ability  to  predict  crack 
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arrest  was  verified  using  the  J-integrai  approach.  The  question  of  a biaxial 
crack  growth  resistance  curve  for  use  as  the  failure  criterion  under  this  type 
of  loading  has  been  explored. 

Six  compLex  structural  panels  have  been  fracture  tested  (two  biaxial 
panels)  and  predictions  made  of  their  residual  strength.  The  predicted  values 
were  in  good  agreement  with  the  experimental  data. 

4.2  PREDICTIVE  PROCEDURE  (STEPS) 

Reference  1 (Section  8.D  detailed  the  steps  involved  in  the  residual 
strength  prediction  method  and  will  not  be  repeated  here  except  in  summary 
form.  1L  should  be  noted  that  the  inclusion  of  biaxial  or  multiaxial  states 
of  stress  have  required  additions  to  the  basic  procedure  (see  e.g. , Steps  3 
and  8).  However,  the  basic  outline  remains  the  same  and  Reference  1 should 
be  consulted  for  additional  details 

Step 

1)  Make  a detailed  finite  element  model, 

2)  Select  crack  lengthts)  for  evaluation. 

3)  Perform  the  analysis  based  on  Dugdale  or  P rand 1 1 - Reu  s s material 
behavior  - stress  state  dependent* 

4)  Determine  stresses  in  the  intact  stiffeners* 

3)  Determine  the  value  of  for  each  applied  stress  to  material  yield 
strength  ratio  (p/Fty)  for  each  crack  length, 

6)  Cross  plot  the  data  of  Step  3 to  form  the  crack  driving  force  curves, 

7)  Determine  the  value  of  normalized  applied  stress  in  the  intact 
stiffeners  at  stiffener  ultimate  strength. 

8)  Obtain  crack  growth  resistance  data  for  the  skin  material  In  uniaxial 
tension  or  biaxial  loading. 

9)  Plot  /Jg  versus  A apj^Y . curve  from  the  data  of  Step  8, 

10)  Overlay  the  data  from  Step  9 on  the  crack  driving  force  curves  of 
Step  b and  determine  the  amount  of  slow  stable  tear  - if  minimal 
(£  0.23  inches)  a stiffener  critical  condition  is  operative  otherwise 
it  becomes  a skin  critical  case  and  both  residual  strength,  slow  tear, 
and  crack  arrest  can  be  determined  from  tangency  conditions. 
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V CONCLUSIONS 


The  following  conclusions  have  been  reached  based  on  the  studies  under- 
taken in  Phase  III  of  this  study.  Following  the  conclusions  a Table  is  given 
in  Section  5.2  which  summarizes  the  accuracy  of  the  prediction  technique  for 
both  Phase  II  and  Phase  III  panels. 

5.1  GENERAL  TRENDS 

The  flexible  fastener  model  developed  in  Reference  1 can  be  used  to 
treat  Load  transfer  in  complex  structural  arrangements  under  complex  loading 
conditions,  Dugdale  material  behavior  can  be  assumed  for  uniaxial ly  loaded 
panels,  P rand tl - Reuss  assumptions  must  be  used  for  biaxial  loading  conditions. 

For  elastic  and  elastic-plastic  analysis*  an  increase  in  square  root  of 
J occurs  under  increasing  biaxial  load,  however  the  stress  transferred  to  the 
intact  stringer  decreases.  Accompanying  this  behavior  is  an  increase  in  crack 
opening  displacement  for  increasing  ratio  of  biaxiality.  This  influence  is 
opposite  of  that  observed  in  unstiffened  sheets.  There  is  a decreasing  trend 
In  plastic  zone  size  with  increasing  biaxial  loads. 

The  procedure  outlined  in  Section  4.2  can  be  used  to  predict  the  residual 
strength  of  typical  aircraft  structure  loaded  in  uniaxial  tension.  The 
extension  of  this  procedure  to  biaxial ly  loaded  panels  has  also  been  explored 
using  certain  limiting  assumptions  with  a good  degree  of  success. 

The  phenomenon  of  crack,  arrest  can  be  predicted  by  using  the  proposed 
residual  strength  prediction  technique.  It  was  determined  that  the  resistance 
curve  of  the  material  must  be  replotted  at  the  new  c rack  length  after  crack 
arrest  occurs. 

5.2  OVERALL  ACCURACY  OF  THE  METHOD 

Table  XV  is  a summary  of  all  structural  panels  tested  during  Phases  11  and 
III  of  this  study.  The  average  difference  between  predicted  and  measured 
residual  strength  for  all  panels  is  seen  to  be  3,1  percent  which  is  well  within 
the  accuracy  expected  of  the  technique  and  well  within  the  materials  data 
scatter  of  ± 5 percent.  As  noted  previously  a prediction  was  not  made  for  the 
biaxial  panel  tested  at  1:2/3  load  ratio  due  to  failure  in  the  attached  doublers. 

Ihe  residual  strength  of  Phase  II  panels  was  predicted  within  10  percent 
for  some  panels,  Ihe  variation  between  predicted  and  actual  failure  stresses 
could  be  due  to  residual  stresses  in  the  panel  accentuated  by  the  chem-mi 1 1 ing 
process.  The  predicted  failure  load  of  the  panel  is  also  dependent  on  the 
assumed  yield  stress  of  the  material.  In  the  analysis  an  average  yield  stress 
obtained  from  coupon  tests  was  used.  The  yield  stress  of  individual  panels 
may  vary.  Therefore,  an  accurate  estimate  of  yield  stress  is  required  for  a 
higher  degree  of  accuracy  in  predicting  residuaL  strength  as  noted  in 
Reference  2. 
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TABLE  XV 


COMPARISON  OF  PREDICTED  AND  MEASURED  RESIDUAL  STRENGTHS 
FOR  PHASE  II  AND  PHASE  111  PANELS 


RESIDUAL  STRENGTH  (ksi) 

PERCENTAGE 

T¥PE  OF  PANEL 

PREDICTED 

MEASURED 

DIFFERENCE 

Phase  11  (Short  Crack) 

•Riveted  Zee  Stiffeners  with  Lands 

40.05 

37.25 

+ 

7.5% 

• Bolted  Zee  Stiffeners  with  Lands 

40.05 

36.39 

+ 

10.1%  . 

•Adhesively  Bonded  Zee  Stiffeners 
with  Lands 

43.8 

42.03 

+ 

4 . 2% 

Phase  II  (Lode  Crack) 

• Rive,ted  Zee  Stiffeners  with  Lands 

27.5 

25.98 

+ 

5 . 85% 

• Bolted  Zee  Stiffeners  with  Lands 

27.0 

26.41 

+ 

2.2% 

•Adhesively  Bonded  Zee  Stiffeners 
with  Lands 

28.8 

28.98 

- 

0.625% 

Phase  III  (UniaxiaL  Tension) 

• Intact  Stringer  - Thin  Skin  (-15) 

39.73 

37.37 

+ 

6.32% 

• Broken  Stringer  - Thin  Skin  (-23) 

34.36 

34.51 

- 

0.44% 

• All  Titanium  Panel 

78.86 

76.45 

3.15% 

• Aluminum  Panel  - Thick  Skin  (-1) 

37.02 

38.78 

- 

4.75% 

Phase  III  (Biaxial  Loading) 

• 1:1/3 

40 

39.87 

+ 

0.33% 

• 1:2/3 

— 

36.98 
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